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FOREWORD 

Results  contained  in  this  report  were  obtained  during  the  course  of 
research  conducted  under  the  contract  N-62269-85-C-0246,  from  the  Naval 
Air  Development  Center  to  Drexel  University.  The  NADC  program  manager  is 
Mr.  Lee  W.  Cause.  The  program  principal  investigator  at  Drexel  University  is 
Professor  A.  S.  D.  Wang.  Contributors  to  this  program  are  Dr.  E.  S.  Reddy,  Mr. 
W.  Binienda  and  Mr.  Y.  Zhong  of  Drexel  University,  and  Professor  F.  Delale  of 
the  City  College  of  New  York.  The  authors  wish  to  thank  Mr.  Cause  for  his 
many  helpful  suggestions  during  the  course  of  this  work. 
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I.  MIXED  MODE  MATRIX  CRACKING 

Introduction- 

Structural  composite  laminates,  notably  those  made  of  unidirectional 
tape  systems,  are  known  to  suffer  extensive  matrix  cracks  before  the  failure 
of  the  load-carrying  fibers.  These  matrix  cracks  occur  because  of  the  inter¬ 
facial  stresses  that  arise  from  local  load- transfer.  For  instance,  load  transfer 
from  a  weakened  fiber  bundle  to  an  adjacent  stronger  fiber  bundle  may  induce 
debonding  between  these  two  fiber  bundles.  Similarly,  load  transfer  from  a 
weakened  lamina  to  an  adjacent  stronger  lamina  may  induce  delamination 
between  these  two  laminae.  Thus,  the  initiation  and  propagation  of  matrix 
cracks  in  laminates  of  multi-directional  plies  usually  follow  either  the 
fiber/matrix  interface  or  the  ply/ply  interface,  or  both. 

Fig  1.1  shows  the  extensive  matrix  cracking  that  ;ccurred  in  a  notched 
graphite-epoxy  [02/902ls  laminate  subjected  to  uniaxial  tension.  Two  major 
types  of  matrix  cracking  occurred  when  viewed  at  the  phenomenlogical  scale. 
Namely,  in-ply  cracking  between  two  adjacent  fiber  bundles  and  delamination 
between  two  adjacent  plies.  The  four  vertical  cracks  were  initiated  first  near 
the  hole  and  propagated  inside  the  0°-plies.  The  driving  force  here  was  clearly 
due  to  the  load-transfer  from  the  fiber  bundle  cut  by  the  hole  to  the  fiber 
bundle  which  is  uncut.  As  the  vertical  cracks  propagated  away  from  the  hole, 
load  transfer  then  took  place  locally  between  the  cracked  0°-plies  and  the 
uncracked  90°-p!ies.  This  in  turn  induced  interply  stresses  along  the  root  of 
the  vertical  cracks;  and  delamination  in  the  0/90  interface  initiated  and 
propagated  with  the  applied  load. 

A  closer  analysis  of  the  cracked  specimen  stress  field  at  each  major 
stage  of  crack  development  would  reveal  a  complex  three-dimensional  stress 
state.  Stress  concentration  is  high  near  the  crack-root  regions;  and  the 
resulting  crack  propagation  is  mixed-mode  (with  both  opening  and  shearing 
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actions!  I]). 

The  basic  mechanisms  of  in-ply  trans-fiber  cracking  and  interply 
delamination  have  recently  been  treated  at  the  phenomenlogical  scale  by  a 
united  energy  method  (2J.  The  energy  release  rate  concept  of  elastic  fracture 
mechanics  was  employed  as  a  criterion  for  crack  propagation.  This  method, 
when  coupled  with  a  three-dimensional  stress  analysis  and  limitted  to  mode- 1 
propagation  conditions,  has  proven  useful  for  predicting  brittle  matrix  cracks 
in  graphite-epoxy  systems.  For  mode- 1  cracking,  it  is  necessary  to  determine 
the  strain  energy  release  rate  G|  as  the  crack-driving  force  and  validate  the 

corresponding  critical  strain  energy  release  rate  G)c  as  a  material  property  12]. 

As  for  matrix  cracks  that  involve  mixed  modes,  such  as  those  illustrated 
in  Fig  l.l,  the  applicabilty  of  the  energy  method  has  not  been  as  firmly 
established.  One  major  difficulty  is  that  the  specimen  geometry  effects  on 
mixed-mode  crack  propagation  are  intrinsically  coupled  with  the  materia) 
effects.  Separation  of  these  effects  requires  a  rigorous  analysis  of  the  cracked 
specimen  at  each  major  stage  of  crack  development. 

There  have  been  several  studies  aimed  at  establishing  criteria  for 
mixed-mode  matrix  cracking  in  unidirectional  laminates.  Most  of  the  studies 
used  graphite-epoxy  systems.  Wilkins,  et.  al.  (3]  and  Ramkumar,  et.al  14]  used 
the  cracked-lap  shear  specimen  loaded  in  uniaxial  tension  to  induce  mixed 
mode- 1 1 /mode- 1  delamination  between  the  cracked  lap-layer  and  the  substrate 
layer.  By  varying  the  thickness  of  the  lap-layer  relative  to  the  substrate  layer, 
the  mixed-mode  ratio  Gjj/Gj  in  the  measured  total  critical  strain  energy 
release  rate  could  be  varied  from  0.35  to  0.45.  They  observed  that  the  total 
strain  energy  release  rate  (Gj*Gjj)c  is  somewhat  greater  than  GIC  obtained 

under  pure  mode-1  conditions.  Bradley  and  Cohen  15)  used  a  cantiliver 


.  ,  v  %  .’j  , 


NADC-87 133-60 


split-beam  specimen  loaded  by  a  pair  of  upward  and  downward  loads  applied  at 
the  tip  of  the  cantlllver.  Variation  of  the  mixed-mode  ratio  6||/G(  was  achieved 
by  changing  the  ratio  of  the  upward  and  downward  loads.  In  this  way,  they 
obtained  a  range  for  the  G||/G|  ratio  from  0  to  about  0.6.  They  observed  that 

the  measured  total  strain  energy  release  rate  (G|*G||)C  increased  uniformly 

with  G|j/G|  in  those  systems  made  of  brittle  matrix  material,  but  decreased 

slightly  with  G||/G(  in  systems  made  of  ductile  matrix  material.  Wang,  et.  al 

[6]  used  a  doubled  side-notched  off-axis  unidirectioal  laminate  coupon  loaded 
In  axial  tension.  By  varying  the  off-axis  angle  from  0°  t  >  90°,  the  mixed-mode 
strain  energy  release  rate  for  Gjj/Gj  ratio  ranging  from  0  to  about  3  could  be 

determined  They  found  that  the  total  strain  energy  release  rate  CGj*G|  ()c  first 

increased  with  G( ,/Gj  monotonlcally  up  to  about  G,,/G,  =  1.5,  and  then  reached 

an  asymptotic  value  which  is  about  three  times  Gjc.  Russell  and  Street  (7]  used 

four  different  test  specimen  configurations  and  obtained  mixed-mode  strain 
energy  release  rates  for  a  wide  range  of  mixed-mode  cracking  conditions.  Their 
results  showed  that  the  measured  strain  energy  release  rate  values  were  test 
specimen  dependent. 

Clearly,  no  firm  agreement  could  be  reached  toward  establishing  a 
general  criterion  for  mixed-mode  matrix  cracking,  despite  these  and  many 
other  efforts.  It  is  believed  that  one  principal  reason  for  the  disagreement  in 
results  is  that  the  effects  of  specimen  geometry  on  mixed-mode  matrix 
cracking  were  not  rigorously  separated  from  material  effects.  Specifically,  the 
mathematical  analysis  for  the  cracked  specimen  was  based  on  either 
one-dimensional  beam  theory  or  a  finite  element  laminated  plate  model.  These 
approximate  methods  lack  the  required  precision  to  treat  highly  three- 
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dimensional  stress  fields  that  often  contain  stress  singularities.  This  is 
especially  true  when  applied  to  mixed-mode  cracking  conditions,  where 
significant  errors  could  result  in  the  computed  6(  and  6|j  energy  release  rates. 

Objectives.  oUtesfcarch 

The  purpose  of  this  study  is  to  overcome  the  above  stated  difficulty  by 
developing  an  analysis  for  mixed-mode  fracture  employing  an  ideal  specimen 
configuration  and  loading  condition.  This  problem  can  be  treated  rigorously  on 
the  basis  of  the  theory  of  elasticity  and  fracture  mechanics.  Concurrently,  a 
finite  element  procedure  is  developed  to  simulate  the  same  problem  and  yield 
accruate  numerical  solutions.  Mathematical  rigor  and  n-./nerical  accuracy  are 
needed  at  the  same  time  in  order  to  ensure  correct  separation  of  the  specimen 
geometry  effects  from  the  material  effects.  Finally,  experiments  are 
conducted  using  specimems  of  similar  configurations  to  generate  mixed-mode 
fracture  data  on  a  wide  range  of  Gj  j/Gj  ratios.  The  corresponding  fracture 

analysis  is  performed  using  the  finite  element  simulation  model.  This  will 
provide  a  final  correlation  between  experiment  and  analysis 

It  should  be  mentioned  that  an  exact  elasticity  solution  for  the  test 
specimen  configurations  cannot  presently  be  obtained,  requiring  use  of  a 
numerical  model  with  established  accuracy. 

Research  Plan 

To  this  end,  the  following  research  plan  involving  four  analysis  steps  as 
illustrated  in  Fig.  1.2  was  implemented. 

Step  A  We  begin  with  a  unidirectional  laminate  of  infinite  domain  Let 
the  laminate  contain  a  through-thickness  kinked  crack  and  subject  to  a  uniform 
fai — field  tensile  stress  o0.  Let  the  base  of  the  crack  be  orientated  normal  to 
the  far-field  tension,  while  the  kink  be  in  the  direction  of  the  fibers  The  angle 
between  the  fiber  and  the  applied  tension  is  denoted  by  8.  Assuming  linear 
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elastic  response  under  the  applied  stress  o0,  a  pair  of  dissimilar  stress 
singularities  develop  at  the  two  crack  tips.  At  some  critical  o0,  the  kink  will 
propagate  in  mixed-mode  along  the  fiber  direction.  Our  plan  is  to  determine  the 

crack  tip  stress  intensity  factors  Kj  and  Kj(  as  well  as  the  strain  energy 

release  rates  Gj  and  Gjj.  By  treating  the  laminate  as  an  elastic,  homogeneous 

and  orthotropic  body  that  contains  the  prescribed  kinked  crack  we  solve  the 
elastostatic  problem  exactly  based  on  the  theory  of  elasticity.  Note  that  the 
determined  mixed-mode  G( j/Gj  ratio  at  the  tip  of  the  kink  will  vary  with  the 

angle  0  as  well  as  the  length  of  the  kink,  a'.  We  are  especially  interested  in  the 
cases  where  a'  is  arbitrarily  small. 

Step  B.  Next,  we  simulate  the  problem  by  a  finite  element  model  Here, 
we  represent  the  infinite  plate  by  a  rectangular  plate  which  contains  the  same 
kinked  crack.  The  length  and  width  of  the  plate  are  so  large  that  the  effects  of 
the  plate  boundaries  will  not  affect  the  stress  state  near  the  kinked  crack.  The 
finite  element  shape  and  mesh  size  selections  must  be  tuned  to  yield  a  result 
as  close  as  possible  with  that  obtained  by  the  exact  elasticity  solutions.  This 
step  establishes  the  degree  of  accuracy  of  the  finite  element  model. 

Step  C.  Having  established  confidence  in  the  accuracy  of  the  finite 
element  model,  we  then  use  the  model  and  simulate  a  finite  width,  off-axis 
unidirectional  tensile  coupon  with  double  side-notches.  We  assume  that  the 
far-field  stress,  o0,  will  initiate  a  kinked  crack  and  propagate  unstably  along 
the  direction  of  the  off-axis  fibers  at  some  critical  value.  The  finite  element 
model  simulates  this  mixed-mode  cracking  and  yields  both  the  kink-tip  stress 
intensity  factors  Kj  and  K(|  and  the  strain  energy  release  rates  G(  and  G(|.  The 

computed  mixed-mode  G|j/Gj  ratio  at  the  tip  of  the  kink  will  vary  not  only  with 
the  off-axis  angle  0  but  also  with  the  depth  of  the  side-notch,  a. 
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Note  that  this  problem  differs  from  the  previous  one  in  that  the  kink-tip 
behavior  in  this  problem  is  affected  significantly  by  the  free  edge  boundary 
and  the  notch  depth  relative  to  the  width  of  the  tensile  coupon. 

Step  D.  The  final  step  involves  experiments  in  which  off-axis  tensile 
coupons  with  double  side-notches  are  tested  the  same  way  as  those  simulated 
in  Step-C  Here,  the  critical  fai — field  stress,  oCr,  at  the  propagation  of  the  kink 
is  recorded.  This  provides  the  necessary  data  for  correlation  with  the  analysis 
performed  in  step-C. 

In  order  to  effect  a  realistic  and  reliable  correlation,  the  off-axis  angle 
0  of  the  test  specimens  is  varied  from  0°  to  90°  in  seven  different  increments 
and  the  initial  depth  of  the  side-notches  is  varied  from  0  1  to  0. 1 75  inch  in  four 
increments  (the  width  of  the  coupon  is  1.0  inch). 

The  analytical  treatment  for  the  problem  as  described  Step-A  and  the 
finite  element  simulation  for  the  same  problem  as  described  in  Step-B  are 
outlined  in  Section  2.  The  detailed  mathematical  derivations  and  solution 
procedures  for  the  analytical  problem  are  included  in  Appendix  A,  while  the 
details  of  the  finite  element  calculation  routine  are  referred  to  an  earlier 
publication.  Section  3  discusses  the  experiment,  the  experimental  results,  the 
finite  element  simulations  and  the  final  experiment/analysis  correlations  as 
described  in  steps  C  and  D  above.  Finally,  a  set  of  concluding  remarks  is 
presented  in  Section  4. 
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Fig  12  Schemalics  showing  the  four  analysis  steps  in  the  present  study 
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II.  MIXED-MODE  FRACTURE  ANALYSIS  MODELS 

In  this  section,  we  shall  consider  the  problem  of  a  kinked  crack  in  an 
infinte  anisotropic  elastic  plate  which  is  subjected  to  a  uniform  tensile  stress 
field.  The  problem  will  be  treated  first  by  an  exact  elasticity  formulation  and 
then  by  a  finite  element  simulation.  The  purpose  of  the  analytical  treatment  is 
to  obtain,  for  the  first  time,  an  exact  elasticity  solution  for  the  posed  problem, 
from  which  the  accuracy  of  the  finite  element  simulation  procedures  can  be 
established. 

Elasticity  Solutions  For  A  Kinked  Crack. 

Consider  the  kinked  crack  shown  in  Fig.  2.1.  The  inifinte  plate  is  a 
representation  of  the  unidirectional  laminate  with  the  fibers  orientated  at  an 
angle  0  from  the  direction  of  the  applied  tension.  The  base  of  the  kinked  crack 
is  situated  normal  to  the  applied  tension,  while  the  kink  itself  is  orientated 
along  the  fiber  direction.  Our  objective  here  is  to  determine  the  mixed-mode 
stress  Intensity  factors  K,  and  K  n  and  the  strain  energy  release  rates  G(  and 

G|j  at  the  tip  of  the  base  as  well  as  of  the  kink. 

The  corresponding  problem  of  a  kinked  crack  in  an  infinte  elastic 
isotropic  plate  has  been  treated  by  Gupta  [8,9]  and  Chatterjee  1 10],  In  both 
cases,  the  singular  integral  solution  scheme  was  employed.  Here,  we  shall  first 
address  the  problem  of  two  separate  cracks  embedded  in  an  orthotropic  plate, 
one  is  the  base  and  the  other  is  the  kink.  Initially,  these  two  separate  cracks 
are  unconnected.  Using  the  crack  surface  displacement  derivatives  as  unknowns 
and  the  principle  of  superposition,  the  problem  is  then  formulated  in  terms  of 
singular  integral  equations  with  Caushy  type  kernels.  The  resulting  system  of 
equations  are  solved  numerically  employing  a  Gaussian  quadrature  and  the 
collocation  method.  The  mode  I  and  mode  II  stress  intensity  factors  are  then 
calculated  with  the  relative  angle,  relative  distance  and  relative  size  of  the 
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two  cracks  as  parameters. 

Next,  by  letting  the  right  tip  of  the  base  and  the  lower-left  tip  of  the 
kink  touch  each  other  it  thus  becomes  one  kinked  crack.  At  this  configuration, 
the  singular  integral  equations  are  still  valid;  but  some  of  the  kernels  become 
singular  at  the  point  of  connection.  By  imposing  the  condition  that  the  surfaces 
at  the  connecting  point  are  now  free  of  stress,  the  problem  reduces  to  that 
posed  in  Fig  2.1.  The  details  of  the  analysis,  including  the  mathematical 
derivations  and  solution  procedures,  are  presented  in  Appendix  A. 

Some  interesting  results  have  been  obtained  for  a  number  of  cases 
Consider  first  the  case  of  two  separate  cracks  where  the  base  crack  and  the 
kink  do  not  touch  each  other.  Their  interaction  in  terms  of  fracture  parameters 
can  be  investigated  by  varying  their  orientation  angle,  relative  length,  relative 
distance  and  the  degree  of  anisotropy  of  the  medium. 

For  example.  Fig  2.2  shows  the  mixed-mode  stress  intensity  factors  at 
the  two  approaching  tips  b  and  c,  as  a  function  of  their  relative  distance,  h,  for 
an  isotropic  plate  in  which  the  two  cracks  are  of  equal  size  (2")  and  v/ith  an 
angle  a  -  30°  (see  insert  in  Fig.  2.2).  The  interaction  between  the  two 
approaching  cracks  is  clearly  illustrated  by  the  dependence  of  the  stress 
intensity  factors  on  the  relative  distance  h.  In  fact,  a  singular  behavior  is 
exhibited  when  the  two  cracks  touch  each  other.  It  is  shown  mathematically 
that  the  crack  tip  singularity  at  this  configuration  is  different  from  being  1/2, 
see  Appendix  A. 

The  corresponding  problem  for  an  orthotropic  plate  (mimicking  a 
graphite-epoxy  UD  laminate)  is  shown  in  Fig  2.3.  The  crack  interaction 
behavior  is  essentially  identical  to  that  shown  in  Fig  2.2 

With  a  fixed  distance  h,  the  interaction  is  expected  to  vary  with  the  kink 
angle  a.  Fig  2.4  shows  the  case  of  an  isotropic  plate  in  which  two  cracks  of 
equal  size  (2‘)  maintain  a  fixed  distance  h  -0  r.  Here,  the  stress  intensity 
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factors  at  the  extreme  crack  tips  (points  a  and  d)  are  plotted  against  the  angle 
a  The  corresponding  plot  for  an  orthotropic  plate  is  sh:  wn  in  Fig.  2.5.  Again, 
the  interaction  behaviors  of  these  two  cases  in  terms  of  fracture  parameters 
are  essentially  similar 

Next,  consider  the  problem  in  which  the  two  cracks  of  equal  size  (2“) 
touch  each  other  at  the  kink  point  and  thus  become  a  kinked  crack.  Fig.  2.6 
shows  the  case  of  an  isotropic  plate,  where  the  mixed-mode  stress  intensity 
factors  at  the  extreme  tips  (points  a  and  d)  are  plotted  against  the  kink  angle 
a.  The  corresponding  plot  for  an  orthotropic  plate  is  shown  in  Fig.  2.7.  The 
general  features  in  these  two  cases  displayed  in  terms  of  the  stress  intensity 
factors  are  also  similar. 

As  was  noted  earlier,  one  problem  of  practical  interest  is  when  the 
length  of  the  kink,  a',  is  shortened  as  to  approach  zero.  In  the  limit,  the 
behavior  of  the  crack  is  still  that  of  a  kinked  crack.  This  problem  is  relevant  to 
the  case  of  a  crack  whicn  is  situated  in  an  orthotropic  plate,  not  aligned  in  any 
one  principal  direction,  but  will  propagate  in  one  of  the  principal  directions 
(eg  the  fiber  direction).  The  mixed-mode  stress  int<  isity  factors  at  the 
initiation  point  of  the  kink  can  be  obtained  by  this  solution  procedure. 

Fig.  2.8  shows,  for  an  isotropic  plate,  the  change  of  the  mixed-mode 
stress  intensity  factors  at  the  kinked  crack  tips  (points  a  and  d)  as  the  length 
of  the  kink  (I2)  is  shortened  to  zero.  Here,  the  base  crack  has  a  length  of  Li  -  2" 
and  the  kink  angle  is  set  at  a  -  30°.  The  corresponding  plot  for  an  orthotropic 
plate  is  shown  in  Fig.  2  9. 

At  this  point,  it  should  be  mentioned  that  the  stress  intensity  factors 
obtained  for  all  cases  relative  to  isotropic  plate  closely  agrees  with  those 
reported  previously,  even  though  our  solution  technique  is  different.  Table  2.1 
shows  a  comparison  between  the  present  results  and  those  obtained  by  Gupta 
(9)  and  Chatterjee  [10]  for  a  kinked  crack  having  the  base  and  the  kink  of  equal 
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size,  and  kink  angles  of  30°,  45°  and  60°.  Thus,  our  results  for  the  isotropic 
plate  can  be  trusted  as  being  exact. 

There  are  no  other  solutions  for  the  orthotropic  plate  to  compare  with 
our  results,  however.  Since  the  isotropic  plate  represent  a  limiting  case  of  the 
orthotropic  plate  in  our  solution,  it  is  believed  that  our  results  for  the 
orthotropic  plate  can  also  be  trusted  as  being  exact. 

Finally,  the  mixed-mode  strain  energy  release  rates  Gj  and  Gt t  at  the 

kinked  crack  tips  will  be  calculated  from  the  solutions  of  the  elastic  stress 
field  which  depends  on  the  solutions  of  the  crack  tip  stress  intensity  factors. 
The  details  of  the  calculation  are  again  presented  in  Appendix  A. 

Fig.  2.10  shows,  for  the  isotropic  case,  the  mixed-mode  strain  energy 
release  rates  G(  andGjj  and  the  total  Gy=Gj*G||  at  the  tip  of  the  kink  (point  d) 

as  a  function  of  the  length  of  the  kink  (L2).  Here,  the  base  crack  length  L 1  =  1  " 
and  the  kink  angle  o  =  30°.  The  corresponding  plot  for  the  orthotropic  case  is 
shown  in  Fig.  2.1 1. 

Fig.  2. 12  is  a  plot  of  the  strain  energy  release  rate\  G(,  G(  (  and  Gy,  for  an 

isotropic  plate,  as  a  function  of  the  kink  angle,  a.  In  this  case,  the  length  of  the 
kink  is  made  small  (L2  =  O.D  compared  to  the  length  of  the  base  (L 1  =  D.  The 
corresponding  plot  for  an  orthotropic  plate  is  shown  in  Fig.  2.13. 

A  more  detailed  parametric  study  on  the  mixed-mode  stress  intensity 
factors  and  the  strain  energy  release  rates  at  the  kink  tip  in  both  isotropic  and 
orthtropic  plates  is  contained  in  ( 1 1 , 1 2). 

Finite  Element  Simulations. 

The  problems,  treated  above  based  on  exact  elasticity  formulation,  will 
now  be  simulated  by  means  of  a  finite  element  model.  The  finite  element 
procedure  employed  has  been  developed  for  a  previous  NADC  study  on  matrix 
cracking  analysis  in  multidirectional  laminates  1 13]  Hence,  the  details  of  this 
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work  will  not  be  presented  here. 

The  simulated  model  is  illustrated  in  Fig.  2.14,  where  a  finite-width 
rectangular  plate  containing  a  kinked  crack  is  considered.  To  simulate  an 

■j 

infinite  plate,  the  length  and  the  width  of  the  plate  are  made  very  large  v 

compared  to  the  size  of  the  kinked  crack.  In  this  particular  simulation,  the  k 

width  W  of  the  plate  was  10  times  the  length  of  the  base  crack;  and  the  length 
to  width  ratio  of  the  plate  was  set  at  3  to  1. 

With  these  provisions,  the  boundary  effects  on  the  finite  element  ‘ 

results  for  the  mixed-mode  stress  intensity  factors  and  the  strain  energy  ; 

release  rates  are  minimized.  Generally,  the  computed  stress  intensity  factors  ! 

are  more  sensitive  to  the  finite  element  shape  and  mesh  size  in  the  crack 
region  than  the  computed  strain  energy  release  rates.  In  particular,  the  total 
strain  energy  release  rate  Gj  ■  G|*G||  agrees  most  closely  with  the  elasticity 

solutions.  Thus,  by  tuning  the  finite  element  shape  and  mesh  near  the  kinked 
crack,  the  accuracy  of  the  finite  element  results  can  be  optimized. 

For  the  case  of  an  isotropic  plate  where  the  kink  is  made  small  compared 
to  the  base  (L2/L1  =  0.1),  a  comparison  is  made  between  the  finite  element 
solutions  and  the  elasticity  solutions  for  the  kink  tip  strain  energy  release 
rates..  Figs.  2.15,  16  and  17  show,  respectively,  the  kink  tip  strain  energy 
release  rates  G(l  G||  and  Gj  as  functions  of  the  kink  angle,  a.  The  corresponding 

comparisons  for  a  kinked  crack  in  an  orthotropic  plate  are  displayed  in  Figs. 

2.18,  19  and  20.  The  close  agreement  between  the  finite  element  results  and 
the  elasticity  results  establishes  the  accuracy  of  the  numerical  procedure. 

The  finite  element  method  will  be  used  to  simulate  the  test  specimens 
used  in  the  experimental  study,  since  an  elasticity  solution  for  those  specimen 
configurations  is  not  available. 
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TABLE  2. 1 

Stress  Intensity  Factors  (normalized  by  o0)  at  Kinked  Crack  Tips 
(Isotropic  Plate  with  Li  •  L2) 


kink  angle 

K/Oo 

Present  Solution  Gupta  (9) 

Chatterjee  (10] 

30° 

Ki(a) 

1.3421 

1.3559 

1  3508 

K2(a) 

0.0328 

00327 

0.0325 

Ki(d) 

1.0949 

1.0873 

1  0830 

K2(d) 

0.6855 

0.6833 

06804 

45° 

Ki(a) 

1.2732 

1.2902 

1.2887 

K2(a) 

0.0217 

0.021 1 

0.0208 

K,(d) 

0.7546 

0.7463 

0.7438 

K2(d) 

0.8450 

0.8405 

0.8377 

60° 

Ki(a) 

1.2082 

1.2221 

1.2194 

K2(a) 

-0.0108 

-0.0109 

-0.01  16 

Ki(d) 

0.394 1 

0.3900 

03822 

K2(d) 

0.8350 

0.8319 

0.8292 
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Fig  2. 1  Geometry  and  reference  frames  for  a  kinked  crack  in  an  infinte 
orthotropic  plate 
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Fig  2  2  Interaction  of  two  cracks  as  a  function  of  their  relative 
distance,  h  Isotropic  case  with  a  -  30° 
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Fig.  23  Interaction  of  two  cracks  as  a  function  of  their  relative 
distance,  h.  Orthotropic  case  with  a  -  30°. 
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Fig  2  4  Interaction  of  two  cracks  as  a  function  of  the  kink  angle,  a. 
Isotropic  case  with  2h/Lj  -0.1. 
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Fig  2.6  Stress  intensity  factors  (normalized  by  o0-/L/2)  at  the  kink  tip 
as  a  function  of  kink  angle,  a.  Isotropic  case  with  U  3  L2 
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Fig  2. 15  Mode- 1  strain  energy  release  rate  (under  eK=  1 )  at  the  kink  tip 


as  a  function  of  kink  angle,  o,  by  elasticity  and  finite  element 
methods.  Isotropic  case  with  L2/L1  *  0  I 
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F  ig  2  1 6  Mode- 1 1  strain  energy  release  rate  (under  e*=  I )  at  the  k  ink 
tip  as  a  function  of  kink  angle,  a,  by  elasticity  and  finite 
element  methods  Isotropic  case  with  L2/L1  -0  1. 
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F ig  2. 1 7  Total  strain  energy  release  rate  (under  ex=  I )  at  the  kink  tip 
as  a  function  of  kink  angle,  a,  by  elasticity  and  finite  element 
methods  Isotropic  case  with  L2/L1  ■0.1. 
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Fig  2.18  Mode- 1  strain  energy  release  rate  (under  e**!)  at  the  kink  tip  as 
a  function  of  kink  angle,  a,  by  elasticity  and  finite  element 
methods  Orthotropic  case  with  L2/L1  "0  1. 
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Fig.  2.19  Mode-ll  strain  energy  release  rate  (under  e*- 1 )  at  the  kink  tip  as 
a  function  of  kink  angle,  a,  by  elasticity  and  finite  element 
methods.  Orthotropic  case  with  L2/L1  -0.1 
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Fig  2  .20  Total  strain  energy  release  rate  (under  e*- 1 )  at  the  kink  tip  as 
a  function  of  kink  angle,  a,  by  elasticity  and  finite  element 
methods.  Orthotropic  case  with  L2/L1  -  0  1 
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III.  EXPERIMENT  AND  RESULTS 
Material  and  Specimens. 

Hercules  AS4-350I-06  graphite-epoxy  unidirectional  tape  was  used  to 
fabricate  8-ply  unidirectioal  (UD)  laminates.  Off-axis  tensile  coupons  9  “  long 
and  1‘  wide  were  machined  from  these  laminates.  End-tabs  1.5"  long  were 
bonded  to  specimens  yielding  a  6-inch  test  section. 

Two  sets  of  coupons  were  prepared  for  test:  Unnotched  and  double-side 
notched.  The  set  of  unnotched  coupons  were  strain-gaged  and  tension  tested  to 
provide  the  basic  unidirectional  ply  elastic  constants.  Coupons  with  the 
following  off-axis  angles  were  tested. 

6  *  0°,  5°,  10°,  15°,  20°,  25°  and  90° 

In  addition,  a  set  of  [*452]s  laminate  coupons  was  also  tested  to 
determine  the  shear  modulus  6LT  of  the  unidirectional  ply. 

The  set  of  double  side-notched  coupons  were  tested  for  mixed-mode 
fracture.  Fig.  3.1  depicts  the  geometry  of  the  double  sic. -notched  specimens. 
Initially,  the  side-notches  were  cut  with  a  8-mil  diamond  saw  to  depth,  a,  as 
shown.  Under  the  applied  tensile  stress  o0,  a  kink  crack,  denoted  by  a’,  is 
expected  to  initiate  from  the  side-notch  tip  and  propagate  in  the  fiber 
direction.  The  propagation  of  the  kink  crack  then  provides  the  desired 
mixed-mode  fracture  condition. 

Since  the  mixed-mode  fracture  behavior  at  the  notch  tip  is  influenced  by 
the  off-axis  angle  0  as  well  as  the  depth  of  the  initial  side-notch,  these 
parameters  were  varied  during  the  experiment  as  follows. 

0  -  0°,  5°,  10°,  15°,  20°,  25°  and  90° 
a  »  O.r,  0.1 25’,  0.15"  and  0.1 75V 

Thus,  a  total  of  28  mixed-mode  fracture  conditions  were  created  during 
this  set  of  tests  Three  to  four  specimen  replicates  were  tested  for  each 
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fracture  condition,  with  the  exception  of  one  side-notch  case  (a=0. 15)  where 
only  one  specimen  was  available  for  some  of  the  off-axis  angles. 

Test  Method  and  Test  Data. 

All  specimens  were  tested  under  room  temperture  condition  in  a 
close-loop  Instron  tester  with  a  load-control  rate  of  4000  Ibs/min.  Tests  of 
0°,  90°  and  [*452ls  unnotched  coupons  were  used  to  determine  the  following 
averaged  ply  properties: 

EL=  21.0  msi  Ey=  1.5  msi  G|_y  =  0.98  msi  t/LT=  0.3 

The  axial  moduli  Ex  of  the  unnotched  off-axis  coupons  were  also 
measured  Their  averaged  values  are  listed  below: 

8 _ 0°  5°  10°  15°  20°  25°  90° 

Ex  (msi)  21.0  18.4  13.5  9.4  6.7  5.0  1.5 

The  critical  load  at  the  onset  of  the  kink  crack  was  recorded  for 
side-notched  coupons.  Post-test  SEM  examination  of  the  f  actured  surface  near 
the  kink  point  was  conducted  to  determine  the  fractured  surface  morphology 
under  500x  to  lOOOx  magnification. 

Fig.  3.2  shows  the  experimental  plot  of  critical  stress  versus  the 
off-axis  angle  0  at  the  onset  of  the  kink  crack  in  those  specimens  having 
side-notches  0.1"  deep,  in  this  plot,  the  average  critical  axial  strains  were 
calculated  using 

(ex)cr  -  (Ox)  cr/Exn  ( I ) 

where  Exn  is  the  axial  modulus  of  the  notched  coupon.  Exn  was  not  measured 
during  the  experiment.  Value  used  in  (1)  was  calculated  by  the  finite  element 
routine.  Table  3.1  lists  the  values  of  Exn  for  all  the  mixed-mode  fracture  test 
cases  conducted  in  this  study. 

Figs.  3.3  and  3.4  show  the  critical  onset  stresses  for  specimens  having 
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notch  depth  a  -  0.125*  and  0.175*,  respectively.  The  case  for  a  -  0.15“  is  not 
shown  because  of  insufficient  numbers  of  test  specimens 

It  is  seen  from  Figs.  3.2  to  3.4  that  the  onset  stress,  (ox)cr,  for  kink 
crack  initiation  and  propagation  decreases  sharply  with  the  off-axis  angle  8, 
but  only  slightly  with  notch  depth.  On  the  other  hand,  the  onset  strain,  (ex)Cr, 
appears  to  be  insensitive  to  the  off-axis  angle  but  sensitive  to  notch  depth. 

Post-test  SEM  examination  of  the  fractured  surfaces  revealed  extensive 
fiber  breaking  near  the  kink.  Fig.  3.5  presents  a  typical  SEM  for  each  off-axis 
specimen  tested  (8  -  0°,  5°,  10°,  15°,  20°,  25°  and  90°).  Fiber  breaks  are 
visible  in  all  cases.  This  phenomenon  is  caused  by  fiber  nesting  near  the  notch  | 

region  and  fiber  bridging  accross  the  cracking  path,  an  indication  of  strong 
adhesion  between  the  matrix  and  fibers. 

Finite  Element  Simulation. 

Using  the  finite  element  routine  discussed  in  Section  2,  the  initiation  ' 

and  propagation  of  the  experimental  mixed-mode  kinked  crack  conditions  can  be 
simulated.  Referring  to  the  off-axis  unidirectional ly  reinforced  composite  I 

coupon  shown  in  Fig.  3.1,  let  the  coupon  be  loaded  by  the  f — field  strain  ex  -  1. 

The  resulting  stress  field  will  be  highly  concentrated  near  each  of  the 
side-notches.  Assume  that  the  notch  tip  stresses  will  cause  a  kink  crack  to  1 

start  from  the  notch  tip  and  that  the  kink  will  propagate  stably  in  the  direction  J 

of  the  fibers.  Of  interest  is  when  the  length  of  the  kink  a’  is  small.  The  j 

mixed-mode  strain  energy  release  rates  G(  and  G||  at  the  kink  tip  are  assumed 

j 

to  control  the  behavior  of  the  small  kink  propagation.  The  values  of  G(  and  G||  - 

■ 

are  calculated  by  the  finite  element  routine  via  a  crack-closure  technique  1 1 3). 

These  are  conveniently  expressed  in  terms  of  the  applied  far-field  strain  in  the  1 

form:  i 

G|-C|(e„)2  G||  -  C, |(e«)2  (2)  J 
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where  C|  and  Cjj  are  numerical  coefficients  from  the  calculation. 

The  values  of  and  the  ratio  for  G(  and  G|(  (or  C(  and  C(|)  can  vary 

dramatically  depending  on  the  off-axis  angle  and  the  sid(  notch  depth.  Fig  3.6 
shows  two  families  of  plots  for  the  coefficients  C|  and  C(|  with  the  side-notch 

depth,  a,  as  a  parameter  and  the  off-axis  angle,  0,  as  a  dependent  variable.  Note 
that  the  kink  is  highly  mixed  in  fracture  modes  for  off-axis  angles  up  to  20°. 
Beyond  25°,  the  fracture  of  the  kink  is  dominanted  by  mode-1.  The  mixed  mode 
ratio,  defined  here  as  Gj j/Gj  and  shown  in  Fig.  3.7,  is  found  to  depend 

principally  on  the  off-axis  angle  8.  It  is  almost  independent  on  the  initial 
side-notch  depth  a. 

Critical  Strain  Energy  Release  Rates  at  Onset  of  Kink. 

The  critical  stress  (oK)cr  at  the  onset  of  the  kink  was  measured 
experimentally.  The  corresponding  critical  strain  (ex)cr  can  be  calculated  by 
mean  of  equation  (1).  These  results  were  shown  earlier  in  Figs.  3.2,  3.3  and  3  4 
for  the  cases  of  a  =  0. 1 “,  0.125*  and  0.175",  respectively.  Thus,  by  using 
equation  (2)  and  the  values  of  C(  and  Cj(  shown  in  Fig.  3.6  we  can  calculate  the 

critical  strain  energy  release  rates  Gj,  Gjj  and  Gj  =  Gj*G||  for  each  of  these 
cases. 

Fig  3.8  shows  the  calculated  critical  strain  energy  release  rates  G|C  and 
Gyc  for  a  »  o  r  as  a  function  of  the  off-axis  angle  8.  If  we  disregard  the 
experimental  scatter  which  is  relatively  small,  the  plot  shows  the  total  strain 
energy  release  rate,  GTc,  to  be  constant  for  all  angles,  dispite  variations  in  the 

mixed-mode  fracture  condition.  The  averaged  Gjc  in  this  case  is  about  1.8 
in-lb/in2. 

Fig.  3.9  shows  the  calculate<N*cain  energy  release  rates  for  the  case  of 
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a  =  0.125".  Again,  the  behavior  of  GTc  is  the  same.  Here,  the  averaged  GTc  is 

about  1.7  in-lb/in2.  The  calculated  Gjc  for  the  case  of  a  =  0.175’  is  almost 

identical  to  that  calculated  for  the  a  =  0.125  case,  see  Fig.  3.10. 

Fig  3.11  is  a  mixed-mode  interaction  diagram  plotted  from  all  the 
mixed-mode  fracture  test  cases,  including  some  from  specimens  of  a=0.15". 
These  data  show  some  degree  of  scatter,  but  the  interaction  diagram  strongly 
suggests  that  a  criterion  based  on  the  total  strain  energy  release  rate  Gjc 

would  be  acceptable.  This  conclusion  is  based  on  the  mixed-mode  fracture  data 
which  consisted  of  nearly  uniformly  distributed  data  points  in  the  Gj|/G|  ratio 

range  from  0  to  2.5.  This  range  had  not  been  achieved  in  any  of  the  referenced 
previous  works. 

Most  available  data  for  graphite-epoxy  composite  critical  strain  energy 
release  rates  are  limitted  to  G|C  Genrally,  the  measured  values  for  G)c  lie  in 

the  range  between  0.7  to  1.5  ln-lb/in2,  depending  on  the  material  system  used. 
In  this  study,  G(c  was  found  to  be  of  the  order  of  1.75  in-lb/in2.  This  seems 

high  compared  to  most  other  accepted  values.  Recall  that  in  all  our  mixed-mode 
fracture  test  cases,  fiber  breakage  or  fiber  bridging  occurred  in  the  wake  of 
cracking,  which  could  account  for  the  higher  measured  values  for  Gjc  and  Gjc. 
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Fig  3.2  Critical  stress  and  strain  at  onset  of  kinked  crack,  a  =  0. 1 
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5  5  (continued)  Photomicrograph  of  fractured  surface  near  Pin):  paint 
A  specimen  of  0  =  10°,  a  -  0  125'  (above),  and  a  specimen  of 
0  ■=  15°  a  =  0  125'. 
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f  k;  3  b  (conhnued)  Photomicrograph  of  fractur eJ  suf  f.K  c  near  Mni  point 
A  specimen  of  0  -  90°,  a  ■=  0  1* 
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Fig.  3.6  Mode-1  and  Mode-ll  strain  energy  release  rate  coefficients  as 
function  of  off-axis  angle  0. 
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Fig  3.7  Mode-11  toMode-1  strain  energy  release  rate  ratio  as  function  of 
off-axis  angle  8. 


Fig  3.8  Critical  strain  energy  release  Rates,  G|C  and  GT<  fora  =  0  1“ 
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IV.  CONCLUDING  REMARKS 

In  this  report,  an  effort  was  made  to  study  the  mixed-mode  fracture 
behavior  of  a  kinked  matrix  crack  in  graphite-epoxy  composite  materials. 
Previous  efforts  in  determining  fracture  quantities,  such  as  the  stress 
intensity  factors  or  the  strain  energy  release  rates,  have  been  based  on  either 
crude  analytical  models  (eg.  simple  beam  theory)  or  finite  element  simulations 
whose  numerical  accuracy  could  not  be  adjudicated. 

(n  this  study,  the  fracture  quantities  K  and  6  were  obtained  both  by 
solving  the  problem  of  an  infinte  orthotropic  plate  containing  a  kinked  crack 
exactly  and  by  a  finite  element  routine.  A  comparison  between  the  exact  and 
the  numerical  solutions  allows  the  finite  element  shapu  and  mesh  size  to  be 
optimized  to  obtain  the  best  results  for  this  class  of  problems. 

Having  done  so,  the  finite  element  routine  is  then  used  to  simulate  a 
similar  class  of  mixed-mode  kink  crack  in  notched  finite  plates,  the  specific 
problem  being  the  off-axis  unidirectional  tensile  coupon  with  double 
side-notches  The  finite  element  analysis  must  be  employed  because  the 
complex  geometry  of  the  test  specimen  would  render  an  exact  solution 
impossible 

in  the  experiment,  we  have  designed  a  specimen  which  is  simple  to  test 
and  is  versatile  in  geometrical  variations  This  provided  some  28  mixed-mode 


fracture  conditions  covering  a  wider  range  of  G(|/6|  ratios  (from  0  to  2.5)  than 

previously  investigated.  Thus,  a  more  definitive  conclusion  can  be  reached 
regarding  mixed-mode  matrix  crack  propagation.  Specifically,  the  test  results 
suggest  that  the  total  strain  energy  release  rate  Gjc  is  a  good  criterion  for 
mixed-mode  matrix  fracture  in  A54-350I  -06  graphite-epoxy  composite. 
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EXACT  SOLUTION  OF  A  KINKED  CRACK  IN  AN  ORTHOTROPIC  PLATE. 

Introduction. 

This  appendix  is  concerned  with  fracture  analysis  of  a  kinked  crack  in  an  infinite  elastic 
orthotropic  plate.  It  is  assumed  that  the  plate  contains  a  through-thickness  crack  of 
initial  length  L-j,  which  makes  an  angle  0  with  one  of  the  principal  direction  of 
orthotropy.  When  the  plate  is  subjected  to  a  far-field  uniform  tensile  stress  normal  to 
the  crack,  the  crack  will  kink  and  propagate  self-similarly  in  the  0  direction.  The  result 
is  a  kinked  crack  propagating  in  mixed  modes,  with  the  degree  of  modal  mixture 
dependent  on  the  angle  0  and  the  ratio  between  the  length  of  the  kink,  L2,  and  the 
length  of  the  main  crack,  L-\ . 

To  determine  the  parameters  relevant  to  mixed-mode  fracture  at  the  tips  of  the  kinked 
crack,  the  problem  is  formulated  in  terms  of  singular  integral  equations  with 
generalized  Cauchy  kernels.  The  resulting  system  of  equations  are  then  cast  in  a  form 
suitable  for  numerical  solution  by  a  Gaussian  quadrature  and  the  collocation  method. 
Analytical  expressions  for  the  stress  intensity  factors,  ki  and  k2,  and  the  strain  energy 
release  rate,  G|  and  G||,  at  the  tip  of  the  kinked  crack  are  also  obtained  in  terms  of 
solution  variables. 

The  problem  is  solved  in  the  following  manner.  First,  we  consider  an  infinite  plate  that 
contains  two  separate  cracks  as  shown  in  Figure  1.  Here,  two  sets  of  reference  frames 
are  used  to  describe  the  two-crack  geometry.  Let  the  applied  far-field  stress  o0  be  in 
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the  y^-direction,  the  base  crack  of  length  L-j  be  on  the  X2-axis,  while  the  kink  crack  of 
length  l_2  be  on  x-|-axis.  The  kink  angle  0  refers  to  the  angle  between  x-|  and  X2-  The 

two  cracks  will  become  one  kinked  crack  if  their  approaching  tips  meet  at  the  common 
origin  of  the  two  reference  frames.  That  is  when  x^  and  xc  become  zero  in  Figure  1 . 


The  two-crack  problem  can  be  represented  by  the  superposition  of  three  individual 
parts,  shown  in  Figure  2.  The  first  part  is  the  plate  without  any  crack  under  the  external 

stress  o0-  The  second  part  is  the  plate  with  only  the  base  crack  of  length  Li  loaded  by 

uniform  normal  stress  acting  on  the  crack  surfaces.  The  crack  surface  normal  stress  is 
equal  and  opposite  the  stress  that  exists  at  the  same  location  in  the  first  part  of  the 
problem.  The  third  part  is  the  plate  with  only  the  inclined  (kink)  crack  of 

A3 
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length  l_2-  Similarly,  on  the  surfaces  of  the  crack  uniform  r.  ma  and  shear  stresses 

are  applied  with  their  respective  magnitude  and  sign  equal  and  opposite  to  that 

0 

existing  at  the  same  location  in  the  first  problem.  The  last  two  parts  of  the  problem  are 
known  as  the  perturbation  problems. 


t  m  m  °o 


mm!  °° 


mini' 

0 


+ 


+ 


perturbation 

problems 


Figure  2.  Method  of  Solution  of  the  Crack  Problem  by  Superposition  of  a  Plate  without 
the  Cracks  Under  Uniaxial  Tension  and  Perturbation  Problems  Under 
Negative  Stress  Imposed  onto  Crack  Surface. 
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Since  the  solution  to  the  first  part  of  the  problem  is  trivial,  w.  are  interested  in  solving 
the  two  perturbation  problems. 

Section  I  treats  the  perturbation  problem  in  which  the  inclined  (kink)  crack  lies  along 
the  principal  axis  (selected  as  the  fiber  direction)  x-j.  The  problem  of  the  base  crack 
situated  along  X2*axis  is  treated  in  Section  II. 

Superposition  of  the  two-part  solution  for  the  two-crack  problem  is  treated  in  Section 
III,  where  a  set  of  singular  integral  equations  are  formulated.  In  Section  IV,  we 
descretize  the  singular  integral  equations  suitable  for  solution  by  a  numerical 
collocation  method.  Expression  for  the  stress  intensity  factors  that  exist  at  the  four  tips 
of  the  two  separate  cracks  are  formulated  in  Section  V. 

The  original  problem  of  a  kinked  crack  is  treated  in  Section  V  by  employing  the  stress 
field  solution  of  the  two-crack  problem  obtained  earlier.  Here,  we  let  the  two 
approaching  tips  to  meet  at  the  common  origin  of  the  xi  -y-j  and  X2-y2  frames.  In  the 
limit  however,  some  of  the  kernels  in  the  singular  integral  equations  become  singular 
themselves.  Conditions  compatible  to  the  new  situation  are  then  formulated,  and  the 
resulting  system  of  singular  integral  equations  is  obtained. 

These  equations  can  be  discretized  in  a  form  suitable  for  solution  by  a  similar 
collocation  method.  Expressions  for  the  kink  tip  stress  intensity  factors  are  thus 
obtained.  This  is  presented  in  Section  VII.  The  definition  of  the  normalized  stress 
intensity  factors  is  given  in  Section  VIII. 
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Finally,  expressions  for  the  mixed-mode  strain  energy  release  rates  are  formulated  in 
Section  IX. 


I.  Formulation  of  the  problem  of  the  crack  along  x-j . 


For  the  inclined  crack,  which  is  situated  along  the  principal  direction  x-j,  the  stress  field 
is  governed  by  the  biharmonic  equation: 

a4F,  a4F,  a4F, 

a22  T  +  (^a22+a66)  2  2+ai1  T =  ° 

ax,  ax,ay,  ay, 

or 


where 


aV  a^F,  aV 

- ^  ♦  P,  — *  P. - 

ax,  ax, ay,  ay, 


(1.2) 


and 


+  a 


22 


66 


aii  "  £ 


LL 


M2 


a22  ~  £ 


TT 


(1.3) 


En,  ^TT'  ^LT  and  VLT  are  the  engineering  elastic  constants  for  the  orthotropic 
material  in  its  principal  frame. 

Using  Fourier  transformation  the  stress  function  is  defined  as: 
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Substituting  (1.4)  into  (1.2),  we  reduce  the  4-th  order  Partial  Differential  Equation 

(PDE)  to  4-th  order  Ordinary  Differential  Equation  (ODE)  with  constant  coefficients  as 

follows:  4  9 

d  *  ,  d  $  4 

P.— *  .-0  (1.5) 

1  «y,  2  <iy2 

The  solution  of  (1 .5)  may  be  written  in  the  form: 

*1(s,y1)  =  A(s)e"“,SVl+  B(s)e'Vyu  C(s) e  D(s)e^yi  (1.6) 

where  coj ,  t02  are  the  roots  of  characteristic  equation 

-  p2  to2  +1  =  0  (1.7) 

where  it  is  assumed  that  Re(co1)  >  0  and  Re(co2)>0. 

By  taking  the  exponetial  behavior  function  for  s  ~  the  stress  function  becomes: 


•M 


A  e 


-“,|s|y, 


Be 


^lslyi]e’isxi 


ds 


Ce“',s|yUDeW 


i  -is*.  -i 

]e  i  ds 


(1.8) 

(1.9) 


where  A,B,C,D  are  functions  of  Fourier  parameter  "s"  to  be  determined  by  the 
boundary  conditions.  Note  that  superscript  of  y-j  "+"  stands  for  positive  y-j  and  for 
negative  y-j. 


It  is  known  that  stresses  are  the  same  on  positive  and  negative  face  of  the  crack.  This 
can  be  stated  in  the  following  mathematical  form: 
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F,(xv0  )  -  F,(xv0“) 
-si-F1(x1,0+)--5i-F1(x1.°-) 


Then,  the  stress  function  becomes: 


f.(xvv!>  *  ^  f  i  a  * B  e^|s,v'i  e‘isx' ds 


(1.10) 

(1.11) 


(1.12) 


^  1 1 K  A  V  B  c2)  eU,|S|v^  *  (A  c3-  B  c,)  •*** ]  e'isx’  ds  (1.13) 


where 


(0  +0) 

c.=  — - — 2.  (1.13a)  c,  =  - 

to  -  to 
1  2 


2to  2co 

- 2_  (i.i3b)  c,= - —  (113c) 

0)  —  to  to  ~  to 

12  12 


Thus,  the  crack  surface  stresses  can  be  obtained  from  the  stress  function: 


»«  =  fs2  [  to.V’W.  EtoV^I  ***'  ds 

¥i  2aJ  1  l  2  J 


(1.14) 


iH=  js2[  aF(Ac,  ♦  Bc2)e“'|s|'C  u=(Ac3-  Bc,)«“zIs|’(  ]  e''“'ds  (1.15) 


oW  =  ±  fs2  ( Ae^W"  ♦  Be  “^1  e'^  ds  (1.16) 

VM  2a  J 

-M 

0H  .ills2  [(Ac.*  Bc2)e“i|s|^+(Ac3-  Bc^Je^Je'^ds  (1.17) 

W  2a  j,  ’  £ 


»W  =  ±  fs!  [  Ae"w,*s|^  +  Be"w^]  e"5"'  ds  (1.18) 

vw  2a  J 

-M 

aH  .  ±1‘  [  (Ac, .  Bc2)e“-M'l  *  (Ac,-  Be,).’**]  e^-ds 
2a  * 


(1.19) 
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A  more  convenient  form  for  the  above  can  be  obtained  in  terms  of  the  crack  surface 
displacement  derivatives  as  unknowns  : 


',(*,)  =  ir  [  u(xvo+)  -  u(xvo')] 


f2(xi)  =  3^-[V(xi.0+)-v(x,>0')] 


for  xc  <  x,  <  xd 


(1.20) 

(1.21) 


where  u  and  v  are  shear  and  normal  crack  surface  displacements  respectively. 

To  relate  the  displacements  and  their  derivatives,  the  material  stress-strain  relation 
(Hooke's  Law)  is  applied  : 


=  a ..a  +  a.„cr 

Yi  11  Yi  12  nH 

(1.22) 

=  a.«cr  +  a.~cj 

ViVi  12  Yt  22  W 

(1.23) 

Upon  substitution  of  the  expressions  for  stresses  (1.14-1.17;,  we  obtain  from  (1.22), 
the  surface  strains: 

-  a12[  Ae‘“»|s|vi  +  Be‘^|s|V|'|}e"s*'ds 

-  ^Js2{a„[(AV  Bc>?d  (Ac,  -  Bc,)co2e  **] 

-  a|2[(Ac,*  Be,)©  °'|S|VL  (Ac,  -  Be,) »  ^K^ds 


SV=  2^  A(a"“' '  a'2>  e'"’WV'  *  B(a"“' '  a'2>  e'”2li|V' 1  e  ds  (124) 


»t.|- tft.n  «  i  «->  n  ii i*i  »'f.  rt  ft  lH  -  rtj  t“'H'rt  «‘>vii 
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J-)  _  1  f ,2, 


2^  Js2[(Ac , -»-  BCjXa,^  -  a,2)  eWl  S  Vl 


2  uJs|y!  -isx, 

+  (Ac3- Bc1)(a11(o2-a12)e  ]e  ds 


(1.25) 


Then,  the  surface  displacement  derivatives  in  (1 .20)  become, 


(1.26) 


Similarly,  by  substituting  (1.14-1.17)  into  (1.23)  we  obtain: 


:y?y, =  h.  K  [A(ai2wi  ■  a22) e  U,|S|y,+  B(ai2w2  -  a22) e  tt2lS|Vl  ] e  ,SXl  ds 


-H  =^}s2[(Ac1+Bc2)(a120)^-a22)e 


®il*|y? 


+  (Ac3  -  Bc1)(a12w2  -  a22)  e  WzlS|yi  ]  e  'S><1  ds 


By  integrating  the  expressions  for  eyiyi  with  respect  to  y-j ,  an  expression  for  the 
normal  displacement  v  can  be  obtained: 


1  Tj,ACai2wt-a22).-“tNy;  B(a,2<Va22)  -«*|s|y;  ,  -isx, 
y  =  — —  s  [ - e  + - e  e  ds 

2ni  -“2|S| 


1  r  2r  (Aci+  Bc2)(a12(o1  -  a22)  «,|s|yj 

__  s  - e 

2nJ  u>1  |s  | 

(Ac3  -  Bc1)(a|2o>2  -  a22)  u  |S|y*  -isx, 
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Now  we  difterentiate  v  with  respect  to  x-j  and  obtain 

1  f  i  S3  r  A(aiiw1  '  a22^  wilsty' 
- e 


av  f  i_s_ 
ax,  2n  J  |s| 


B(a)2u)2-a22)  -ujlslv;  -•**, 
- e  e  ds 


U)-, 


(1.27) 


av  _  -i 
ax,  2n  J 


i  s3  (Ac,*  Bc2)(a,2uj,  -  a22)  «,|s|y; 

-  -  <3 

|S  I  1  «, 


(Ac]  -  Bci)(a120)2  -  a22)  o  |s|yj  -is*, 
- e  1  e  ds 


Substituting  (1 .24)  and  (1  25)  into  (1  26),  and  (1.27-1.28)  into  (1 .21),  we  obtain: 


(1.28) 


1  f  2  i  -«,|s|y,  *»,l*|y,  j 

ft(Xi)=hm  _  s  |[Ae  -  (Ac,+  Bc2)  e  ](anu>,  -  a,2)  * 

y,-*0  L u  ^ 

-[Be  •  (ACj -  Be,) e  IV|  ](at|o>| •  a,2))  e  ^  ds  (1.29) 


yx,)=l.m  AjiiJ 
y,-»0  J,  I  I 


Ae  +  (Ac,+Bc2)e 


u> 


„  *«2|s|yj  „  «»|s|yj 

Be  *(Ac3-Bc,)e  2  1 

+ - (ai2w2 '  a22)f 


-|S*i  ^ 

e  ds 


(1.30) 


or 


r’(Xl)“  2n  J 


s2((A-Ac1-Bc2)(a11u)^a)2) 


4  (B  -  Ac3  *  Bc))(anco2-a,2)]e  'S>S  ds 


(1.31) 


A  1 1 
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, ,  1  M  s3  r  (A  +  Aci  +  ®c2)(ai2a,i  *  a2 2) 

*2  *1' —  2n  J  Is  I  1  T 


(B  +  Ac3- Bc,)(ai2o>2- a22)  -isx, 
- ]  e  as 


(1  -32) 


Upon  substitution  in  (1.31-1.32)  for  c-j ,  03,  C3  from  (1.13  a-c),  the  following  is 


obtained: 


fi(Xl>=  Iri  1  2s2a11(Aa>i+  Boj^o)^  u>2)e  ^  ds  (1  -33) 


-If  is3  ta.+  a)^  -isx. 

w- :d2^(A*B>±^e  ds 


(1.34) 


CD  CD 
1  2 


The  inverse  Fourier  transformation  of  (1 .33)  and  (1 .34)  gives: 


1  f  is,i 

A cd  +  Bco_  =  2  I  ^  M  e 

*  2  2a,,(“,*  “2)  s  i 

A  *  B  -  -£L - -  [  tj  t>)  dt, 

is  2a22(u> 


W1+U)2  >  v 


with  A  and  B  taking  the  forms: 


A  =  — 


Vt  ♦  4-^V-r  'iw*' 

s  2a,1(u)1  -  co2)  ^  i  S  2a22(u),  -  u>2)  ^ 


-1  1  f  ■'*»,  ls(  r  ist, 

B--~"  2-  2.  Iw  «.-H  -V  2,  [w  «. 

s  2a11(a>,  - u>2)  ^  IS  2a22(wr (U2)  ^ 


(1.35) 


(136) 


Equations  (1.35)  and  (1.36)  can  be  substituted  in  (1.14),  (1.16),  (1.18)  such  that  the 


VVV.V  ,V,V.  ./  s“.-/  V  V.  V.\  V.V.V.V  V  ».■  *.  •  -  -.  V  v 


#«%'«■  9  ^  1^* 
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crack  surface  stresses  are  expressed  in  terms  of  f-j ,  fa  After  integrating  with  respect  to 
"s"  and  applying: 

2  2  ^22 
i  2  a1t 

these  stresses  of  the  perturbed  problem  are  obtained  as: 


1 


f1(t1)^yl*f2(t1)(t,-x1)0i, 


\{ 


^  2tl(w2- W2)an  ^ 


2  2..  .2 

“iVl  *  (t,  xi> 


•)  2 


(1.37 
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II.  Formulation  of  the  problem  of  the  crack  along  X2- 

For  the  horizontal  crack,  which  is  situated  in  the  X272  ^rame  anc*  is  making  the  angle 
0  with  the  direction  of  fibers,  the  governing  field  equation  is  a  general  biharmonic 
equation  for  anisotropic  materials  in  the  form: 
a4F, 


a4F„ 


aV 


22  4 

ax2 


-  2b 


26  3 

ax2ay2 


+  (2b12+b66) 


22'  2bt6 


ax2ay2 


4  4 

a  F,  a% 

— \  +  t>„-T=  0 

ax2ay2  ay2 


where 


bn  =  ancos40  +(2a12+  a66)sin28  cos20  ♦  a22sin48 

4  2  2  4 

b22  =  a, ^in  0  +  (2a12+  a66)sin  0  cos  0  +  a22cos  0 

bl2  “  a!2  *  (a„*  a22  "  2ai2  -  a66)Sir,2e  CO!;20 

b66  *  a66  +  <ai  I  *  a22 ' 2ai2  ‘  a66)Sin20  COs20 
2  2  1 

bi6  =  la22sin  0  *  ancos  9  +  -^(2a)2+ a66)cob  28  ]  sin20 
b26 =  [az2c0S^9  *  ansin20  -  -^<2a12+a66)cos28]sin20 


Let's  define: 


Pi-* 


2b 


26  . 


22 


P2  ” 


^b12+  b66 


22 


p3  =  " 


2b 


16 


22 


b11 

P<-  e11; 

VZ2 


and  after  substitution  of  (2.2)  into  (2.1 ),  the  4-th  order  PDE  is  reduced  as  follows: 


aV, 


ax. 


a4F, 


ax2ay2 


P 


*'F2 

2ax2ay2 


P 


J^2_ 

3  3 

ax2ay2 


+  P<- 


*2 

*A 


Assume  the  solution  for  the  stress  function  F2  in  the  following  form: 

A14 


(2.1) 


(2.2) 


(2.3) 


I 


$ 

2 

V 

V 

V 

V 

V 


I 


A 

A 


,v 


a 

1 


V\J 
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F2(x2.y2)-^12ck(s)eW®"Vds 


then,  the  tollowing  characteristic  equation  is  obtained  by  substituting  (2.4)  into  (2.3): 


-  ip3r3  -  p2r2  +  i  pt  r  +  1  -  0 


The  roots  of  (2.5)  are  of  the  general  form: 


r,  =  a  +  i  b  ;  r2=c+idj 
r4  =  -  a  +  i  b ;  r3  =  -  c  +  i  d; 


where  a,c  >  0  is  assumed. 

Thus,  for  y2  >  0,  the  stress  function  F2  can  be  written  as: 


f2(x2'  Yi)  “  2n  {  C1  e 


(a^y^  (MidjyjS  -  ix<s 

+  C3  e  ]  e  ds 


+  [C2e 


(-c-wdjy^  -ix^ 


+  e 


]  e  ds 


Define  new  constants  as  follow: 


.  fc,  rs< 

|  C2  if  s  > 


(2.8)  B(s) 


J  C3  if  s<( 
"  1  C4  if  s  > 


then  (2.7)  becomes: 


1  r  -•|s|y*+iby*s  -c|s|y’+idy^  -bcj* 
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Similarly,  for  y2  <  0,  we  have: 


(-a-HbJy^  (<-Hd)yzs  -ix-s 

C2  e  +  C4  e  ]  e  ds 


f  (a-Hb)y2s  (c-wdjy^  -iiy  1 

+  I  [  C1  e  +  C3  e  ]e  ds  J- 
o 


(a-HbJy^  (c-HdJy^  -ix^ 


(2.11) 


We  then  define  the  following  constants: 


C(s)  = 


=  r  c2 if  s< 

\  C.  if  s  > 


(2.12)  D(s)  = 


_JC<  if  s< 

l  C3  jf  S> 


(2.13) 


and  (2.11)  becomes: 


r,v  -v  1  frr^S|y2  +  iby^  n  clslV2  +  idV2!t  ‘ix2s  H 

f2(x2,  y2)  =  J  l C  e  +De  Je  ds 


(2.14) 


Since  the  stresses  are  the  same  on  the  positive  and  negative  face  of  the  crack  we  can 
state: 


F2(x2,  0  )  =  F2(x2,  O’) 


(2.15) 


^f2(V0+)»^f2(x2,0-) 


(2.16) 


Satisfying  the  above  conditions,  the  stress  functions  (2.10)  and  (2.14) 
become  respectively: 


= -k  ii  a  e 


•|s|y2-t-iby^s  -c|s|y*+idy^  -ix^ 


+  B  e 


]  e  ‘  ds  (2.17) 
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/c  A  O  DN  elsly^  +  idy^s  -ix^ 

+  (S3A  +  SAB)e  *]e  ds 


(2.18) 


where 


c  _  |s|(c  +  a)»ls(d-b).  r  - 2 a |s| 

1  IsKc-al  +  isfd-b'i  '  K> iy >  b3 "  kiTr— T7TT 


|s|  (c  -  a)  +  i  s  (d  -  b) 

S2=  |s|(c-a)  Jis(d-b)'  (2,19b) 


|s| (c -  a)  +  is(d-b) ;  (2  19c) 

-isjcc^iscd-b) 

4  Js|(c-a) +  is(d-b)  v  1 


The  second  derivatives  of  the  stress  function  give  us  the  crack  surface 
stresses: 


-i}|A(-a|s|*tos)V*|s|y^lbv^ 


x2xz  2n 


c/  i  i  \2  -clsly2+idy^  -ix^ 

+  B(-c|s|  +  ids)  e  1  *]e  ds 


J-)  _  ±  f  2  r  /c  a  *  o  D.  ^ly^iby^ 

2n  JS  i  + 


y2y2 


c|s|y,  +  idyls  -ix,s 

+  (S3A  +  S^B)e  12  V2j  e  ^  ds 


(2.20) 


+  (S3A  +  S4B)(c|s|  +  ids)2e  1  j  e  ***  ds  (2.21 ) 


4* 

(+)  ‘1  f  2 . .  **|s|yi  +  iby^  -e|s|y*+idy^s  -ix_s 

°vtyt"arJs  ,A®  *B®  i«  <*  (2.22) 


(2.23) 
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^y2  =  ^:|  ist(S,A  H.  S2B)(a|s|  *  ibs) e**s|y2+,bVzS 

c|s|yl+ idyls  -ix-s 

+  (S3A  +  S„B)(c|s|  +  ids)e  2  r]e  ^ds  (2.25) 

Similarly,  as  in  Section  I,  the  discontinuity  at  the  crack  is  expressed  as  a 
set  of  crack  surface  derivatives: 


f3(  x2)  *  TT I  U(X2'  0+)  *  u(x2'  °")1  (2-26) 

2  for  x  <  x  <  x 

f4(  X2)  =  I  Y(X2'  0+)  -  V(X2*  °*)  )  (2.27) 

The  stress-strain  relations  (Hooke’s  Law)  in  the  X£-y2  frame  are: 

c  =  b.«<j  +b.«o'  +b.tT 

*2*2  *2*2  ^  ^2^2  ^  X^Z 

c  =  b^o  +  b-.cr  +  b.,x 

V&Z  ^  x2X2  ^2  Yz^Z  ^  *2^ 

Upon  substitution  of  the  expressions  for  stresses  (2.20-2.25)  into  Hooke's 
Law,  the  surface  strains  are  obtained.  Then,  integrating  the  £^2  with 
respect  to  y2>  an  expression  for  the  normal  displacements  is  obtained, 
which  is  differentiated  with  respect  to  X2  and  substituted  into  (2.27).  The 
expression  for  ex2X 2  030  be  directly  substituted  it  .a  (2.26).  Finally  the 


following  expressions  are  obtained: 
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|bii|A[(-alsl+ibs)2-  S,(a|s|  +  ibs)2  -  S3(c|s|  ♦  ids)2] 
b[(-c|s|  +  ids)2  -  S2(a|s|  +  ibs)2  -  S^(c|s|  +ids)2 1 


The  inverse  Fourier  transformation  of  (2.28-2.29)  gives: 

S 


AP,  +  BP2  -  |f3(ye'Sl*dt2 


AP3+  BP,  =  { 


from  which  we  find  that: 
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where 


P,  =  s2Q,;  P2  =  s2Q2;  P2 

=  sQ3 

;  p< 

Q*  =  2b,,[a(a  +  c)-  ia  (b-d)] 

for 

s>0 

Q"j  =  2b, , [a  (a  +  c)  +  ia  (b  -  d)] 

for 

s<0 

Q*  =  2b, ,[c  (a  +  c)  +  ia  (d  -  d)] 

for 

s  >  0 

Q2  =  2b, ,[c  (a  +  c)  -  ia  (b  -  d)] 

for 

s  <  0 

(2.33) 


(2.34) 


,2  2 


2a[(db  -  d  -  c  -  ac)  +  i(cb  +  ad)]  u 

U3  "  .2  2.  .  2  2  b22 

(a  +  b  )(c  +  d  ) 


.2  2 


2a[(d  +  c  -  db  +  ac)  +  i(cb  +  ad)]  u 

3  -  2  ~2~  2  .2  b22 

(a  +  b  )(c  ♦  d  ) 


for  s  >  0 


for  s  <  0 


(2.35) 


2  .  2 


-  2c[(a  +  b  -  db  +  ac)  +  i(cb  +  ad)] . 

Q< - -  ,  2  h2v  2  V - bZZ  f0r  S  »  0 

(a  +  b  )(c  +  d  ) 


-  2c[(db-a2-b2-ac)  +  i(cb  +  ad)] . 

“  *  ,  2  2  -  bZZ  f°r  S<° 

(a  +  b  )(c  +  d  ) 


(2.36) 


Substituting  the  above  for  A  and  B  into  (2.20),  (2.22)  and  (2.24)  the  stress 
field  is  obtained  : 


1  'ff  «,  W  ■  ■ 
2a  I  ^  y2(a+ib)  + 


’  R2  W  (  R3  W  ■  R<  ‘3^) 
+i(t2*x2)  +  V2(c+id)+i(V  x2) 


,  r7  W  -  Re  W' 

y2(a-ib)  - 1(1^-  x2)  +  y2(c-id)-i(t2-x  ) 


(2.37) 
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i  r  r  ^9^^) '^10^2)  (  r,iw-f 

:  2a  J  I  y2(a+ib)  +i(L,-  x2)  +  y2(c+id)  + 

**1  ■ 


l^(^)  ~^i2^(^) 

(c+id)  +i(^-  x2) 


.  RnW-W  rhW-W1  , 
y2(a-ib)  -  i^-  x2)  +  y2(c-id)  -  ifl-  x2)  J  dt2 


(2.38) 


T  _  _1_  [  T  ^17^3^)  *^18^(^)  R19^( 

^2  2  ft  J  I  y2(a+ib)  +1(1^-  x2)  *  y2(c+ 

*1  ■ 


y2(c+id)  +i(l^-  x2) 


^  -  *yA> .  w  -  yA)!  Jt 

*  y2(a*'t>)  -  i(tj*  X2)  +  y2(c-id)  -  i(t>-  x2)  I  dt* 


(2.39) 


where  R's  are  coefficients  in  the  form: 


R,--J(a.ib)2; 

°5 


R, - f  (a  *ib)z ; 


I-Oi,  2 

R3  = - (c  +  id)  ; 

°5 


2 

R4  «  —  (c  +  id)  ; 

°5 

°4  2 

Rs  =  “7  (■  a  +  ib)z; 

Qs 

iQ2  2 

R6  =  '~('a  +  ib) ; 
^5 


IVJ.  J 

r7 - -(-c  +  id) ; 

#"VT 


U3  2 

Rg-— (-c  ♦  id)  ; 
Qc 


r9=~; 

°5 


R,3  =  +  ' 

Q5 


Vv ‘vv»-vv V v  wV\  ' 


Vvv;.i  viTtVviiJ 
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iQ 


R1? - ^(a  +  ib); 

CL 


Q2 

R18  =  4(a  +  ib); 
Qs 


Q 


R19  =  — (c  +  id); 


Q 


5 

iQ! 


R20  =  — (c  -  id); 
CL 


*Q2 

ri.  =  -7; 

Q5 


iQ 


R,s'< 


R16  =  -  — ; 

q; 


iQ 


R2i=“7(-a  +  ,b); 
Qc 


°2 

R^2  —  +  (*  a  +  ib) , 

Q5 


QJ 

R23  =  —  C-C  +  icO; 
°5 

ic>; 

— 7(-«  +  W); 
° 


Note  that  the  stresses  in  the  first  problem  (1.37-1.39)  can  be  written  in  the  form 
analogous  to  those  in  the  second  problem  (2.37-2.38): 


j_7f 

Vi  =  2nJ|  y^+i^-Xj)  +  y,<V  i(V  x,) 

"*"5  ^l^l)  *  "^6  ^"7  ^2^l)  *  L, 

+ - + - I  at, 

y,<v  i(Vxi)  y,w2-'(Vxi)  J 


(2.40) 


A  22 


TV  v  >,'VJ.".V.V.V.V.V„V.VV. 
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j_f,rT.w-TiA(ti)lT..w-v.( 

2nJ  I  y^i+i^-x,)  y,w2  ♦  i(t,  -  x, 

■ 


T»W- W2E1)  ,  T15r2(tl)~T16fl(tl)1 

yi»i-Kt,-Xi)  y1w2  -  i(t,  -  x,) 


(2.41 ) 


Vi  =  2a  J  I  y,a)t+  i(t,-  x,)  +  y,o)2  ♦  i(t,  -  x,) 


T2,W- v^) ,  v^) 

vW’^-x,)  ^iw2~ 


M]dt) 

-  Kt,  -  x,)  | 


(2.42) 


2(a)2 -u)2)an 


2(0)2 -0)2)  a, 


T9‘T13“T12e=TU 


2(0)2  _  qj2 
v  1  2 


*  T,0  -  T14  *  ~  .  2  2v  * 

20^0)2-  0)2)  an 


2(0)2 -0)2)  a, , 


2u)2(0)2-0)2)an 


2(0)2 -0)2)  a, , 


2  (0)2 -0)2)  a, 


"^18  “  ^22  =  ^19  =  "**23 


2(0)2 -0)2)  a  t1 


2  (o)‘  -  o)2)  a 


A  23 


jwuuuuuuul1  juuwlw  wuuuw.wulv  'jj  'J  ^  ww!*m  w 
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III.  Generation  of  Integral  equations. 


WW.v.v.v 


By  superposition  of  stresses  from  the  two  problems,  we  obtain  solutions  for  a  plate  with 
two  cracks.  Thus,  we  proceed  with  the  coordinate  transformation  in  order  to  express 
the  total  stresses  in  one  common  coordinate  system.  If  the  total  stresses  are  expressed 
in  xi  -y-|  coordinate  system,  the  transformation  is  given  by: 


(3.1) 


x2  =  x,  cosB  -  y,  sin0 
y2  =  x,  sin0  +  y1  cos0 

and  the  stresses  of  the  second  problem,  when  expressed  in  x-|-y-|  system  are  given 
by: 

°)CX  +  CTv»v»  CT>cx»  ’  CTv«v. 

(3.2) 


-fallaa .  !s^c0S2e  * ,  « 

xixi  2  2  ^2 

dy  «  --  VzV2  -  ****  ^  cos 20  -  T  sin20 

MM  £  t  ^>2 

t'  =  -  VzV2  sin28  ♦  t  cos20 

*1yi  2  *2^ 


(3.3) 

(3.4) 


and  the  total  stresses  are  given  by: 


T 

/r  e  #T 

+  a ' 

*i*i 

O  •  Q 

*i*i  *1*1 

T  =CJ 

4.  re  ' 

Vi  Vi  Vi  Vi 

VlVl 

tT  -r 

4.  *r  * 

Vi  Vi 

Vi 

(3.4a) 


Similarly,  if  the  total  stresses  are  expressed  in  X2-y2  coordii  Je,  the  transformation  is 
given  by: 

x,  =  x2cos0  +  y2sin0 

(3.5) 


x,  =  x2  cos0  +  y2  sin0 

y,  =-x2sin0  +  y2cos0 
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and  the  transformed  stresses  of  the  first  problem,  when  expressed  in  X2*y2  are  given 

— -2— v?  *  ^4^ cos2e  ‘  Tx,yisin20  <3-6> 

gv2v2 *  ~~~ cos 2 8  ♦  TXiViSin20  (3.7) 


<j  -a 
*i*i  v,v, 


sin28  ♦  t  cos26 


Thus,  the  total  stresses  in  X2-y2  system  are: 


a  ■>  a  *  a ' 

*2*2  *2*2  *Z*2 


<J  «  ff  +  (J ' 

y2y2  y2y2  v2y2 


T  =  t  +  T' 

*2V2  *2V2  V2 


(3.8a) 


Consequently,  the  transformed  stresses  of  the  second  problem,  when 
expressed  in  x-j-yi  frame,  (3.2)  to  (3.4)  have  the  form: 


a’  =-^~ 


H,  W  +  H2  W 


(x^inO 

+  y^osSXa  +  ib) 

*■  -xtcos0 

+  y1$in0) 

H3  r&) + 

"SA) 

(XjSin© 

+  y^osOXc  +  id) 

+  i^-XjCosB 

-►  y1sin0) 

H5  W  * 

H6  W 

(XjSinO 

+  y1cos8)(a  -  ib) 

-  i^  -  xtcos0 

+  y^inB) 

h7  w  ♦ 

H.W 

1 

(x^inB  +  y,cos0)(c  -  id)  -  i(t^  -  XjCos©  +  y^inQ) 


f 


where 


and 
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H,  =  j[  R,  +  R9  +  (R,  -  R9)cos26  ♦  2R1?sin29  ] 

H2  = '  jl  R2  +  Rio+  (R2  ‘  RiO)cos20  +  2R18sin20  ] 
H3  =  -  j(  R4  +  R12+  (r4  *  R12)cos28  +  2R2Osin20  ] 

H4  =  y[  R3  +  R1t  +  (R3-  Rn)cos28  +  2R1Qsin29  ] 

h5  =  j(  R5  +  Rt3+(R5-  R,3)cos20  +  2R21sin20  ] 

H6  -  -  J [  R6  +  R14+  (R6  -  r14)cos20  4.  2R22sir^0  ] 

H7  -  -  7 1  R8  +  R16+  (R8  -  R16)cos28  ♦  2R24sin20  ] 
1 

Hg  ~~2  [  R?  +  Ri5+  (r?  ‘  Ri5)c°$20  +  2R23sin20  ] 


1 


ll 


H5  W  ♦  HlA(y 


yiyi  2n  I  I  (x,sina  *  y,cos8)(a  ♦  ib)  ♦  i(tj  •  x,cos8  ♦  y.sine) 

*1 


HnW  ♦  Ht2f<(y 


(x^ine  +  y,cos8)(c  +  id)  +  i(^  -  XjCose  +  y^inO) 


HofjCtj)  *  H„%) 


(x,sin0  ♦  ytcos0)(a  -  to)  -  -  x^osO  +  y,sin0) 


H,5yy  *  h,,^) 


(x,sin0  ♦  y,cos0Xc  -  id)  -  i(tj  -  x,cos8  *  y,sin0) 


] 


<*2 


(3.10) 


a 


i 


Yj 

i 

& 


V? 

» 

5 


:  V 


I 

i«v 


V  J 

S 


\  ^ 

:‘-0! 

*  .  ^ 

i 

l 

v.v 

'.K* 

VN 
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where 

Hg »  j  [  R1  +  Rg  -  (R1  -  Rg)cos28  -  2R17sin29  ] 
Hio“  ‘  7l  R2  +  Rio'  (R2 '  Rio)cos28  -  2R18sin28  ] 
Hi  i = '  7 1  R4  +  R12"  (R4  *  r12  )cos20  -  2R20sin26  ] 
h,2=  \[  R3  +  Rn  - (R3  -  Rn)cos28  -  2Rtgsin28  ] 

HJ3  =  -y[  R5  +  Ri3“(R5*  r13 )c°s28  -  2Rzlsin20  ] 
H„-  -  j[  R6  ♦  Rm-  (R6  -  R14)cos20  -  2R22v :n20  ] 
Hjj  =  ■  Rg  +  R|5—  (Rj  *  R]j)cos28  -  2R2<Jsin29  ] 

i 

H16=  — [  Ry  +  R^+  (R7  *  R15)cos20  +  2R23sin20  ] 

and 


_±7T  _ "ilW  _ 

\v\  ~  2a  J  I  (XjSine  +  y^osBXa  +  ib)  +  -  x^osB  +  y1sin0) 

HuW  +  H2oW 

4.  - 

(x^inQ  +  ytcos0)(c  +  id)  +  -  XjCOsS  +  y^inS) 

^2ivy +  ^22  V2) 

+ - 

(x,sm8  +  y,cos8)(a  -  ib)  -  -  x,cos8  +  y,sin8) 


(3.11) 


_ H23W  _ 1  & 

(x,sin0  +  y^oseXc  -  id)  -  i^  -  x(cos8  +  y ,sin8)  I  ^ 
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where 

Hi7“^R9‘Ri)sin20  +  2R17cos20] 
Hu“^(R2*Rio)sin20  ‘  2Rtgcos28 1 
H19“^(R4-Ri2)sin2e  -  2R2Ocos20] 

H2o“^i(R,rR3)sin20  +  2R19cos20] 
H2,  =|[(R13-R5)sin20  +  2R21cos20] 
H22  =  ^(R6-RH)Sin20  -  2R22COS20] 
H23“^(R8-R16)Sin20  -  2R24COS20] 
H24  ■  -^{(R^  *  R?)sin28  +  2R23cos28  ] 


Similarly,  the  transformed  stresses  of  the  first  problem,  when  expressed  in  X2-y2 
frame,  (3.6)  to  (3.8)  have  the  form: 


-J-?  r  °iwto2w 

2n  J  I  (-x2sin0  +y2cos8)u1  +  i(t,  -  x2cos0  -  y2sin0) 

G3  W  +  G4  W 

♦  - 

(-x2sin0  +y2cos0)u>2  +  i(t,  -  x2cos0  -  y2sin@) 

Gs -  GeW 

4.  _____ - - — - - 

(-x2sin0  +y2cos0)u>,  -  i(t,  -  x2cos0  -  y2sin0) 

<hW  *  G«  W 

+  "  — ~ 
(-x2sin0  +y2cos0)w2  -  i(t,  -  x2cos0  -  /2sin0) 


(3.12) 
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G 

G. 

i 

i,  -  j[ T,  *  T,  ♦  (T|  -  T,)cos2e  -  2T17sin28  ] 

!  -  •  |lT2  *  T.o  +  (T2  •  T,o)cos20  -  2T„sir>281 

G 

3  -  •  \  [ ' \  *  T12  ♦  (T„  ■ -  T12)cos28  -  2T20sin28 1 

G 

«  -  y[  T3  *  T„  *  (T,  -  T1()cos20  -  2T„sin28  ] 

G 

5  =  —  [  T5  +  T13  +  (T5  -  T13)cos20  -  2T21sin20  ] 

G 

6---jiTe,*Ti<’(V  T„)cos28  -  2T2jsin28  ] 

G 

,  -  41  T»  *  T'‘ * (T«  ■  vcos2e  ■  2T2<sin29 1 

Q. 

g  =  —  [T7  +  T1S  +  (T?  -  T15)cos20  -  2T23sin20  ] 

*< 
i  r 

r  Q»w  *Qu>w 

1™ 

H 

1  (-x2sin0  +y2cos0)a)1  +  i(t1  -  x2cos0  -  y2sin0) 

G„f,(t,)  *Q„W 

(-x2sin0  +y2cos0)co2  ♦  i(t,  -  x2cos0  -  y2sin0) 

0., W  ♦<W.) 

(-x2sin0  +y2cos0)u),  -  i(t,  -  x2cos0  -  y2sin0) 

Q»W 

(-x2sin8  +y2cos0)w2  -  i(t,  -  x2cos9  -  y2sin0) 
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where 


and 


G9  =  2[  W  <V Tg)cos20  +  2T17sin20 ] 
Gjo  =  *  +  TJ0  -  (T2  -  Tlo)cos20  +2T1gsin20] 

Gu  *  *  "2  [  T4  +  T^2  -  (T^-  T12)cos20  +2T2Osin20  ] 

Gi2=  ^  T3  +  Tii  - (T3  -  T1t)cos20  +2T19sin20  ] 
G13=  2"  [  ^"5  +  "^13  ~  O5  *  T13)cos20  +2T21sin20  ] 

gh  =  4[T6  +  Ti4  ' (Te * tm)cos20  +  2T_  sin20  ] 
G15  = "  "2"[  ^8  +  *^"i6  ~  ("^8  *  T,6)cos20  ♦  2T24sin20  ] 
G|6=  ~2  ["^*7  +  ^"is  “ (^"7  *  "Tis)cos20  +2T23sin20  ] 


ll 


qitW  +  gisW 


V2  2ti  J  I  (-x2sin0  +y2cos0)m1  +  i(t,  -  x2cos0  -  y2sin0) 

\  ■ 

Q19f,(t1)  *G20f2(t,) 

+  - 

(-x2sin0  +y2cos0)u2  +  i(t,  -  x2cos0  -  y2sin0) 

G2lf,(t,)  ♦VA) 

+  - 

(-x2sin0  +y2cos0)w1  -  i(tt  -  x2cos0  -  y2sin0) 

G23ft^l)  +  G24^l) 

+  - 

(-x2sin0  +y2cos0)tD2  -  i(t,  -  x2cos0  -  y2sin0) 


] 


dt. 


where 


Gi7  =  lKVT9)sin20  +  2T17cos20] 
Qi.--j!(T,o-T2)sln2e  -  2T1(cos29  I 

G.S-il(Tl2-T<)Sin2e  -  2T20COS2e) 
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G20-^I(T3-T,,)sm2e  .  2T]9cos28  ] 

Q2,-^[(T5-T13)sin20  ♦  2T2,cos28] 

Ga-jKWinM  •  2T22cos28] 

G23  =  ^ t(Ti3 "  T2)sin28  -  2T„cos20] 

G24  =  7KVT,5>sln20  *  2T2Jcos28] 

So,  the  total  stress  field  for  the  infinite  plate  with  two  cracks  is  expressed  by  the  stress 
components  either  expressed  in  x-j  -y-j  frame: 

’JTViW-TiW  Wi)-VA) 


H,  W  ♦ 

HzW 

(x^inB  ■ 

*■  y1cos0)(a  +  ib)  ■ 

f  i^-x^cosB 

+  y^sinB) 

H3  W  + 

H«W 

(XjSine 

♦  ytcos8)(c  +  id) 

+  i^-XjCose 

+  ytsinB) 

H5  W  * 

(x^ine 

+  y1cos0)(a  -  ib) 

-  -  xtcos8 

+  ytsin8) 

H7  W  + 

H.W 

(x^inS 

+  ytcos0Xc  -  id) 

-  i(t>  -  x,cos8 

+  y,sin8) 

V 
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[Wi)-TioW ,  TnW- Wi) 

yiO>1+  i(trx,)  +  y^2  ♦  i(t,  -  x,) 

Wi)-W,)  TMV-TMVl^ 

+ - + - I  dt. 

yi«1  -  '(t,  • x,)  y,iu2-  i(t,-xt)  J 


Hg  W  + 

(x^inS  i 

■  y^osSXa  +  ib) 

+  i^-XfCosB  h 

-  y1sin0) 

"uW  ♦ 

HizW 

(x^inS  -i 

►  y,cos0)(c  +  id) 

+  i(^  *  x^osO  •< 

^sinS) 

^13^)  + 

"MV 

(x^ine 

+  yJcos0)(a  -  to) 

-  i(t>  -  x^osB 

+  ytsin0) 

H«W  + 

h  MV 

i 

(x^inS 

+  y^osOXc  -  id) 

-  I(^>x1cos0 

+  y^inS) 

v  v  -■  v 
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Sv,  =  2a  J  I  y^+  i(t,-x,)  +  y,o>2  +  i(t,  -  x,) 

\  ii 


^2/1^1)  ~  “^22^1 )  +  1~23^l)  ' 

y^ri(trx,)  y1a)2-i(trxi) 


H,7f&>  +  H1.W 


(x^ine  i 

•  y^osBXa  +  ib) 

+  i^-XfCosB 

+  y^inB) 

^19  ^(*2)  + 

H2oW 

(x^inB  1 

i-y^osBXc  +  id) 

+  i(t^  -  x^osO 

+  y^ine) 

H21W  + 

H22W 

(XjSinB 

+  y1cos0)(a-  ib) 

-  i(L  -  x^osO 

♦  y^inS) 

H23W  + 

"zaW 

(x^ine 

*y,cos0Xc  -  id) 

-  i^  -  x^osB 

♦  y^mO) 

(3-17) 


or  the  X2-y2  frame: 


A  33 
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(-x2sine  +y2cos8)u>1  +  i(t,  -  x2cos9  -  y2cin9) 

_ G3  +  G<  tfo) _ 

(-x2sin8  +y2cos8)w2  +  i(tj  -  x2cos0  -  y2sin0) 

_ Ggf^,)  ♦  Q6f2(tt) _ 

(-x2sin8  +y2cos0)co1  -  i(t,  -  x2cos0  -  y2sin0) 


_ G7f1(t1)  +  QgW _ 

(*x2sin0  +y2cos0)co2  -  i(t,  -  x2cos0  *  y2sin0) 
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t  _1_  r  I"  +  RnWRiA(^) 

Vz“  2n  J  I  y2(a+ib)  ♦i(t2-x2)  +  y2(c+id)  ♦Kt2-x2) 

■ 

+  y2(a-ib)  -  x2)  +  y2(c‘id)  *  i(t2‘  x2)  J  ^ 

jjr  Q,r,(t,)  .g1oW 

+  2n  J  J  (-x2sin0  +y2cos0)u)1  +  i(t,  -  x2cos0  -  y2sin0) 

G,iW+GiaW  . 

+  -  —  ■  -  —  1  » ■ 

(*x2sin0  +y2cos0)oj2  +  i(t,  -  x2cos0  -  y2sin0) 
Gia^i)  +G1/2(t1) 

+ - 

(*x2sin0  +y2cos0)a)1  -  i(t,  -  x2cos0  -  y2sin0) 

g«W  +  Q16r2(t1) 

+ - 

(*x2sin0  +y2cos0)u)2  -  i(t,  *  x2cos0  -  y2sin0) 


i 


(3.19) 


I 
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T  1  f  [  (  ^19^^'^20^3( 

Sv* "  2fi  J  I  y2(a+ib)  +1^-  x2)  +  y2(c+id)  -h^- : 

^  w 


*«  r 

bl[ 

*e  L 


tW&WA)  ¥A)-Wl  d, 

’  y2(a-ib)  -  K*2-  x2)  +  y2(c-id)  -  i^-  x2)  J  ? 

_ Qi7W  *  ^18^1) _ 

(-x2sin6  +y2cos0)w,  +  i(tt  -  x2cos0  -  y2sin0) 

_ ^19  *  ^20*2^1) _ 

(-x2sin0  +y2cos0)a>2  +  i(t,  -  x2cos0  - y2sin0) 

_ Q2i^(^i)  ^22^1) _ 

(-x2sin0  +y2cos0)(D,  -  i(t,  -  x2cos0  -  y2sin0) 

_ Q23^i^i)  +  ^24^1) _ 

(-x2sin0  +y2cos0)w2  -  i(t,  -  x2cos0  -  y2sin0) 


(3.20) 


In  the  stresses  (3.15-3.17)  or  (3.18-3.20),  there  are  four  unknown  functions  f-|,  ^  f3-  U 
which  must  be  determined  from  the  appropriate  boundary  conditions  for  this  problem. 
Since  the  plate  is  under  uniaxial  tension  a0  ,  for  the  crack  lying  on  X2  (xa  <  X2  <  x^) 


we  have: 


(3.21) 


v° 


(3.22) 


And,  for  the  crack  lying  on  x-|  (xc<xi  <  x^)  we  have: 


v,--°ocos2e 


^  -aosin0cos0 


(3.23) 


(3.24) 
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applying  these  four  boundaries  and  normalizing  the  variables  x  and  t  by 
substituting  the  following: 


..  xz(\  ■  x»)  xb  +  x. 

2  4  2 

(3.25) 

s2(x„  •  *.)  *b  *  *„ 

X2  “  2  +  2 

(3.26) 

dt2*  2  dx2  -1<t2,s2<1 

(3.27) 

and 

.  X,(Xd-Xe)  Kd»Xc 
^  “  2  4  2 

(3.28) 

Sl(Xd'Xc)  Xd+Xc 

xi  2  +  2 

(3.29) 

xj  ■  xe 

dt^-^dt,  -1  <t1/Si  <  1 

(3.30) 

a  set  of  Cauchy  type  singular  integral  equations  is  obtained  as  follows: 


C12  +  1  K13f3(T2)dt2  +  I  K14f<t(T2)dT2  =  -°0COS2Q 

,  T,  -  S,  1  .  * 


(3.31 ) 


-1  '1  -1 


C2,  |  -i^-dT,  +  J  K23f3(T2)dT2  *  j  K2,f4(T2)dT2  -  -ffosin8cos0  (3.32) 

«  "  Si  .  * 


-1  '1  1 


CJ3  1  7^"dT2  4  CM  \  TT6'!  *  IKl>WdT1  4  jK32f2(T|)dT1  -  -o  (333> 


-1  2  2 


-1  '2  2 


^43  |  dt2  +  ^44  j  -  dt2  +  }^41  ^1^Tl)  dt1  +  j^42  dt1  “  0  (3  34) 

'  T«  *  So  ,  t>‘  So  ,  , 


-1  '2  2 


-1  '2  2 
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which  must  be  solved  for  the  unknown  functions  f-j  .f2»f3  and  U-  For  single-valuedness 
of  solution  it  is  necessary  to  use  also  additional  conditions  which  in  this  case  are  going 
to  be  equations  of  the  following  form: 

i  i 

|f1(x)dT  =0  (3.35)  }f2(T)dT  =  0  (3  36) 


J  f3(  x )  dx  =0  (3.37)  }^(T)dx  -  0  (3.38) 

-1  *1 

The  constants  C's  in  (3.31  -3.34)  are  real  numbers  and  are  calculated  according  to  the 
following  expressions: 


C12  = 


T14  +T11  "T10'Tt5  . 

c  »R"- 

R10  +  R14  ‘  R15  . 

2ai 

°34 

2d 

T17  *  T20  *  T21  +  T24  . 

c  -R"- 

R20  ”  R21  +  R24  . 

2ai 

C43~ 

2rci 

R0  -  R12  -  R13  +  R,6 

r  R19* 

R18  +  R22  “  R23  . 

2ni 

U44* 

2tii 

The  Cauchy  type  kernels  K;;  are  defined  as  follow: 


[s,(xd  -  x)  *  xd  ♦  x  ][sin8(a  *  ib)  -  icosSJ  +  i[t  (x.  -  xj  *  x  +  xj 


[Si(xd  -  xc)  ♦  xd  ♦  xc][sin8(c  ♦  id)-  icos8)  +  i[;2(xb  -  xj  +  xb  +  xj 


[s,(xd  -  Xc)  +  xd  +  Xc][sin8(a  -  ib)  +  icosB]  -  i[x2(xb  -  xj  *  xb  +  xj 


[st(xd  -  xe)  +  xd  ♦  Xc][sin8(c  -  id)  +  icosB]  -  i[x2(xb  -  xj  ♦  xb  +  xj 
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[s,(xd -  Xc)  +  Xd  ♦  xj[sin0(a  +  b)  -  icos0]  ♦  i[x2(xb  -  xj  +  xb  ♦  xj 


[s,(xd  -  xc)  +  xd  +  xe][sin0(c  +  id)-  icos0]  ♦  i[x2(xb  -  xj  +  xb  +  xj 


[s,(xd  -  Xc)  +  xd  +  Xc][sin0(a  -  b)  +  icosOJ  -  i[x  (x  -  xj  +  xb  +  x  ] 


[s1(xd  -  xc)  +  xd  +  xc][sin0(c  -  id)  +  icos0]  -  i[x2(xb  -  x  )  +  xb  +  xj 


Is, (Xd  -  Xc)  +  Xd  +  XcKSin8(a  +  ib)  -  iC0S8l  +  '[Tz(Xb  -  Xa)  +  Xb  +  XJ 


Is, (xd  -  xc)  ♦  Xd  ♦  Xc][sin0(c  ♦  id)-  icos0]  +  i[x2(xb  -  xj  +  xb  +  xj 


ls,(xd  -  xe)  +  xd  +  xe][sin0(a  -  b)  +  icos0]  -  i[x2(xb  -  xj  ♦  xb  +  xj 


[S,(Xd  *  Xc)  +  Xd  +  Xc]tSin8(C  '  id)  +  iC0S8l  -  '[X2(Xb  -  Xa)  +  Xb  +  XJ 
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[s,(xd  -  xe)  +  xd  ♦  xe][sin0(a  +  ib)  -  icosB]  ♦  i[x2(xb  -  xj  ♦  xb  ♦  xj 


[s,(xd  -  xc)  +  xd  ♦  xc][sin8(c  ♦  id)  -  ico$0]  ♦  i[x2(xb  -  xj  «■  xb  ♦  xj 


[s,(xd  -  xc)  ♦  xd  +  xc][sin0(a  -  ib)  ♦  icos0]  -  i[x2(xb  -  xj  +  xb  ♦  xj 


(s,(xd  -  xc)  +  xd  +  xj[sin0(c  -  id)  ♦  icos0]  -  i(x2(xb  -  x 


m] 


-  [  S2(Xb  '  XJ  +  Xb  +  XJSin0  "l  +  iCOS0l  +  iITl(Xd  *  Xc)  +  Xd  +  XJ 


[  s2(xb  ■  \)  +  \  +  X Jsin8  W2  ♦  icosB]  ♦  i[x,(xd  -  xc)  +  xd  ♦  xj 


-  [  s2(xb  -  x J  +  xb  +  xj[sin0  ui  -  icos0]  -  i[x,(xd  -  xc)  +  xd  ♦  xj 


-  [  Sz(Xb  *  \)  *  Xb  +  XJ(S'n0U,2  -  iCOS0]  -  '(Tl(Xd  '  Xc)  +  Xd  +  Xcl 


iwuj 
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K„- 


VXc 


32  2 n  I  -  [  s2(xb  -  x J  ♦  xb  ♦  x  JsinB  u>,  ♦  icosB]  ♦  i[t,(xd  -  xc)  ♦  xd  +  xj 


( s2(xb  -  x J  +  xb  ♦  x  JsinB  u>2  +  icosB]  ♦  i[r1(xd  -  xe)  +  xd  +  xj 


*  I  s2(xb  -  xj  ♦  xb  ♦  x  JsinB  «oi  -  icosB]  -  |x,(xd  -  xe)  ♦  xd  ♦  xj 


- 1  s2(xb  -  \)  ♦  xb  +  x JsinB w2  -  icosB]  -  i[r,(xd  -  xc)  ♦  xd  +  xj 


= 


VXc 


41  2n  I  -  [  s2(xb  -  xj  +  xb  +  x JsinB  w,  +  icosB]  ♦  i[x,(xd  -  xc)  +  xd  +  xj 


[  s2(xb  *  XJ  +  xb  +  xJIsin0(O2  +  icos8l  +  '[Ti(xd  -  xc)  +  xd  +  XJ 


*  I  s2(xb  ‘  XJ  +  xb  +  xj[sin0a>i  -  icosB]  -  i[r,(xd  -  xc)  ♦  xd  ♦  xj 


- 1  s2(xb  -  xj  ♦  xb  +  x  JsinB  a)2  -  icosB]  -  i[t,(xd  -  xc)  ♦  xd  +  xj 


K„,= 


Xd’Xc 


42  2ft  I  -  [  s2(xb  -  xj  +  xb  +  x JsinB  u,  +  icosB]  +  i[r,(xd  -  xe)  4  xd  4  xj 


-  ( s2(xb  -  \)  *  xb  +  xJsjn0  w2  +  icos0l  +  'tTi(xd '  xc)  +  xd  +  XJ 


[  s2(xb  -  xj  4  xb  4  x JsinB  -  icosB]  -  i[x,(xd  -  xc)  4  xd  4  xj 


-  ( s2(xb  -  xj  4  xb  4  xjsin8w2  -  icosB]  -  |x,(xd  -  xc)  4  xd  4  xj 
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IV.  Method  of  solution  -  set  of  Cauchy  type  integral  equations. 


The  set  of  equations  (3.31-3.34)  are  known  as  singular  Cauchy  type  integral 
equations.  They  can  be  solved  by  using  Lobatto-Tchebyshev  collocation  method.  In 
this  method,  we  replace  each  equation  by  n  linear  algebraic  ,quations  of  the  following 
structure: 

F  (t^)  Wk 

C„„2 -C— —  *  2  F<\) wk  -  P(*j>  (4.1 ) 

1  k=1  Tk-Xj  k=t  "p  1  1 

where  a,p  =  1 ,2,3,4  and 


W)  =  W"  ~  2(n- 1) 
Xi  C0s  2n-  2 


k=1,2,...,n 

(4.2) 

i 

w,-n.1  r-2-3 . rvi 

(4.3) 

-“.I 

■0 

_  ■ 

J-  1.2 . n-1 

(4.4) 

*V 

v 

,s' 

V 

.V 

Also,  equations  (3.35-3.38)  are  replaced  by: 


n 

Iwv° 


Hence,  the  set  of  integral  equations  together  with  conditions  are  reduced 
mathematically  to  set  of  4n  algebraic  equations  with  4n  unknow;, in  the  form: 

[A!  (F)  -  {  P  )  (4.6) 

where  [  A  ]  is  4n  x  4n  matrix  of  coefficients  and  {  P  }  is  loading  function  vector.  So  the 
solution  of  the  problem  is: 

{  F  }  *  |  A  1 {  P  }  (47) 

A  42 


(4.7) 
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V.  Stress  intensity  factors. 


Stress  intensity  factors  at  the  tips  of  the  crack  are  defined  as  follow: 


kt(xj  « lim 

s/2(Xm  -  h)  V2<X2'  0) 

(5.1) 

k2(xj  =  lim 

sl 2(X*-X2)  ^Vi(x2<0; 

(5.2) 

kl(Xb)  °  lim  V 

*2 -**6 

^2(X2  -  Xb)  <Jy2y2(X2'0) 

(5.3) 

k2(xb)  =  lim 

n/2(X2  -  Xb)  Tx^y2(X2/0) 

(5.4) 

k,(xe)  -  lim  N 

/2(xe-x,)  °y1v,(xl'0) 

(5.5) 

k2(xc)=lim  > 

x)~*\ 

/2(xc - x,)  V|(*,,0) 

(5.6) 

k1(xd)  =  lim  , 

s/2 (X1  ‘  Xd)  ffyiyi(Xl20) 

(5-7) 

k2(xd)  =  lim  , 

s/2(x1-x<J)  x^y^x^O) 

(5.8) 

where  the  stresses  are  the  total  stresses  (the  superscript  T  is  omitted). 

It  is  also  possible  to  express  the  stress  intensity  factors  in  t  rms  of  a  special  function 
evaluated  at  the  ends  of  the  crack  ie.  F^(1)  or  F^(-1)  (fi=1,2,3,4)  by  applying  the 
following  transformation  to  the  singular  part  of  stress  equations. 

It  is  known,  [1]  and  [2],  that  the  singular  behavior  around  the  end  points  of  the  crack 
can  be  expressed  by: 
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1  [  M.dt  -  -  F(a)(b  -  a)*  * — -(2  -  *  F(b)(b  -  a)» 

n  *  t_z  sinrr(3 


sintux 


'  (2  •  b)*  +  0(z)  (5.9) 


where  a,p  are  order  of  singularity  at  a,b  respectively;  z  is  a  complex  number;  0(z)  is  a 

higer  order  error  term  which  can  be  neglected.  Function  F(t)  is  assumed  to  be 
bounded  because  singular  behavior  of  f(t)  is  completely  determined  by  the 
fundamental  function  w(t).  For  this  class  of  problems  in  asticity  the  fundamental 
function  is  of  the  form: 


w(t)=  (b-t)“(t-a)fl 


(5.9a) 


Thus,  we  can  write: 


f  (t)  S  F_( t )  (b -  t)”(t - a)n 


(5.9b) 


If  the  singularity  at  the  crack  tips  is  -0.5  and  a=-1 ,  b=1  equation  (5.9a)  can  be  written 


f  (t)  = 


F(t) 


(5.9c) 


Using  the  above  equations,  it  is  assumed  that,  the  order  of  singularity  at: 


(5.9d) 


It  will  be  shown  in  the  case  of  a  skewed  crack  configuration,  that  the  order  of 


singularity  at  the  common  point  (x5=xc=0)  is  not  equal  to  -0.5,  while  at  xa  and  x^  it 
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remains  -0.5. 


Applying  the  transformation  (3.26), (3.29)  into  equations  (5.1 -5.8)  we  obtain  the  stress 


intensity  factor  ki  (xa)  as: 


Vx. 


0)  ■ 


Vx. 


sz~*  -1 


■  -~2_” xi +  s2)  WV')  = 


(5.10) 


=  lim  .  JVA  a  (S2,0) 


S2-»  -1 


All  the  singular  terms  in  equations  (3.31-3.34)  are  transformed  using  (5.9).  For  the 
case  of  a=-1  and  b=1 ,  the  first  two  terms  of  (3.31 )  become: 


nC« { 1 777 dT2 •  ■ nC»  ( • 1 -hr 's2  * 

-1  2  S2  X  Sin(-Tr) 


1  e 


\  ^  sin(-y) 


(5.11) 


i  f  dT*  =  f  F<(‘1  > 7f  -77  1 )' 

-i  t2  s2  X  v*  sint-y) 


F<(1)7sin(VS2'1,Z 


(5.12) 


Finally  substituting  the  sum  of  (5.11)  and  (5.12)  for  Cy2y2(S2.0)  into  (5.10)  we  have: 


W-1\/^  [C«F3(’1)  *  C»F71)] 


(5.13) 


VI.  Skew  Crack  Configuration  Xfc=xc=0. 

Introduction. 

By  setting  Xb=xc=0  in  equations  (3.15-3.20)  we  have  the  case  of  a  skewed  crack  in  an 
infinite  plate.  The  boundary  conditions  (3.21-3.24)  remain  the  same  and  equations 
(3.15-3.30)  become: 

A  46 


v.v 

-1 

v.'Ji 

H 
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*2-0 

(6.1) 

x2  =  (1  -s2)y 

-1  <  T2<s2  <  1 

(6.2) 

dt2  =  *y  dx2 

(6.3) 

tt=(1  +Ti)y 

(6.4) 

xi  =  (1  +si)y 

-1  <TvSt<  1 

(6-5) 

x, 

dt,  =  ydTi 

(6-6) 

Normalizing  the  integral  equations  with  (6.1 -6.6)  we  obtain  again  a  set  of  integral 
equations  as  (3.31-3.34),  however  the  kernels  become  singular  when  t^  and  Xj 
approach  the  common  point  (xb=xc=0).  Thus,  a  different  kind  of  collocation  method 
has  to  be  employed.  Also  (3.35-3.38)  no  longer  holds,  so  another  set  of  constraints 
must  be  applied. 

If  one  attempts  to  use  the  Lobatto-Tchebyshev  collocation  method  as  a  rough 
approximation,  then  the  following  conditions  can  be  used: 


F1(-1)=  0 

(6.7) 

f2(-D  =  o 

(6.8) 

o 

u 

/*> 

- - - 

r> 

Li- 

(6.9) 

F<(1)  =  0 

(6.10) 

since  we  have  to  assume  that  the  function  f(t)  is  approximated  by  F(t),  for  example: 


iwwwu«BMfv»wywi  rx  ■>  vm  xv^y.  ^  - 
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! 

i 

C* 

V' 

i 


FjTO 

f,C3  s 


but  as  it  was  mention  above,  in  the  skew  crack  case  at  the  common  point  (%=xc=0) 
the  order  of  singularity  is  not  -0.5,  the  proper  fundamental  function  is  actually: 


F3(t) 

(j(t)  S  -*L=(t-  1)- 


where  F  is  new  bounded  function  and  a  is  not  -0.5.  Considering  the  above  we  have: 


F,(t)  -  F^/TTft-l)*  =>  F3(1)  -  Fj(1)  0  =■  0 


ave: 

similar  arguments  (6.8-6.10)  can  also  be  proved. 

It  is  better,  however,  to  use  the  Gauss-Jacobi  or  Lobatto-Jvtoobi  collocation  method, 
because  the  general  order  of  singularity  p  can  be  applied.  In  this  case,  we  can  use  the 

additional  constraint,  that  the  displacements  at  the  common  point  must  be  the  same, 
either  for  the  main  or  for  the  kink  crack: 


k. 


xd  {  w  <*,  -  x.  j  [  W  SinB  '  w  cose  ]  <*2 


-1 


*d  I  W dx1  -  -  X.  I  [  VT2>  COS0  *  W  SinB  ]  dx2 


-1 


(6.11) 


(6.12) 


We  also  need  two  additional  equations  and  to  determine  the  order  of 


singularity  p  at  the  "wedge  appex"  ie.  common  point  of  both  cracks  (0,0). 
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Let's  define  new  variables  as  follow: 


P2  =  t2'1  *2 <  p2 <  ° 

dp2  =  dt2 

(6.13) 

and 

w,  *  [cos 8  +  isin0  (a  +  ib)] 

(s,  *1) 

(6.14) 

w2  =  [cos0  +  is  in  8  (c  +  id)] 

<■?> 

(S,  *1) 

(6.15) 

w3  *  [cos 8  +  isin8  (a  -  ib) 

I(S,  *1) 

(6.16) 

*  [cos8  +  isinB  (c  -  id) 

1(5,  *1) 

(6.17) 

Noting  structural  similarity  between  K 

13.K14.K23  and 

K24  also 

substituting 

(6.13-6.17)  in  (3.31)  and  in  (3.32)  we  obtain: 

r  i  f  W  H  i 

nC,2n  J  t  .s  dT»  2 
-1  T1  b1 

* 

‘H* 

■ 

1  f  W 
Ri  P2'W1 

dp2 

i  r  f3<p2) 


*H  ~  f 

"fl  J  o  .  w 

-2  P2  2 


dP2  *  h,4  I 


1  f  ^P2) 


P2-W3 


dp2 


-2  P2 


1  f  W 


P2'W1 


dp2 


(6.18) 


'H|2  J  n..w.  ?2  *Hm«J  p2.w3  ?1 


P2-W2 


1 


1  f  f4(P2)  .1  2a 

- dp  I  -  -cr  cos  0 
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and 


nC2ilfW  i|  ,H  if  w_d 
!I«!t.-S,  '  2  I  '7«dp2-w, 


-1  '1  *1 


,  ,  1  f  f3<P2)  ^  (J  1  r  f3(P2> 

Hi^Jr^rdp2+H2iaJ 


P2-W2  ‘  "n±P2-*2 


dp. 


..If  f3(P2)  „  L1  1  r  '4^ 

+  H23—  - dP2  ■  Hi«  r  - dP; 

23fljp2-w<  ttjp*-w1 


1  r  W 


„  1  ?  W  „  „  If  W  , 

H20n  J  rt  .w  dp2  +  H22  xi  J  -  dp2 

■2 


P2  '  W2 


1  f  f^(p,)  | 

H,.,— - dp,  I  *  -<Jftsin8cos8 

24  *  p.-w.  H2  I  0 


P2-w3 


P2‘W4 


] 


Performing  similar  operations  on  (3.33)  we  obtain: 


Ksi  2a  I  (-—  )(1  -s2)(sin8<D1  +  icosB)  +  i(t,  +  1) 
I  X  . 


_ 22U _ 

(-  — )(1  -  s2)(sin0  w2  +  icos8)  +  i(T,  +  1) 


'13 


(-—Jo  -s2)(sin0u)t  -  ic°s0)  -  1(1,  +  “I) 

xd 

Q, 


*15 


(•  ^rJc1  •  s2Xsin0  u2  - icos0)  ‘  '(T1  +  ^ 

"J 


'j  ■  y  f  ijj 


(6.19) 
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or 


+ 


+ 


-iG, 

(t,*1)  -  il 

■ 

X  I  X 
o-  \P 

(1  -SjXsinew,  +  icos8) 

•iGn 

(T1+1)-  i 

(-—  )(1  -$2)(sin6w  ♦  icos6) 

Xd  2 

iG13 

(t,  ♦  1)  ♦  il 

(4) 

*d 

(1  -  s2)(sin0  to,  -  icosB) 

iGis 

(t,  ♦  1)  ♦  i| 

[-  — )(1  -  s2)(sin0  w2  -  icosB) 

Again  defining  new  variables: 

0<p,<2 


2,  =  [cos8-  isinB  w,] 

(4) 

*d 

/  K\ 

z2  =  [cos0-  isin0cj2] 

'(v1) 

z3  =  (cos8  +  i  sin0  w,] 

(•£) 

*d 

'(Sj-1) 

zA  -  (cos8  +  i  sin0  u>2] 

(-*) 

'(S2- -t) 

Substituting  (6.20-6.24)  into  (3.33)  and  (3.34)  we  obtain: 
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It  should  be  noted  that  for  s-|-*  -1  when  -1  <  s-j  <  1  Re(w1,w2,W3,W4)  >  0  and  for 

S2~»  1  when  -1  <  S2  <  1  Re(z-j  ,22,23,24)  <  0.  So  applying  (5.9)  to  (6.18)  and  noting 
that  terms  with  {3  and  f4  are  singular  only  for  p2=Wj=0  or  *2=1  and  Si  =-1  we 
have: 


nCi2  f-  F2(-1)  2*f  cotnp,  (s,  ♦  I)”’  ♦  F2(1)  2'  cotnoc,  (1  -s/' 


F )  2*  — —  w"*  +  0(p)  1 
J  sinrux^  I 

F3(1)  2°*— ! —  w**  *  0(p)  1 
simto^  *  I 

FJ(1)2"1— w‘*.0<p)  1 
sintux^  I 

0.  1  «,  I 

F3(1)2  *— L_w  *+0(p) 
sinoo^  1 

f),  1  «,  1 

F4(  1 )  2  — — — w,  +  0(p) 
sinoo^  I 

0,  1  1 

F^(1)  2  *— w2‘+0(p)  I 
sinru^  I 

F'0)2'— —  wj‘*0(p) 1 
Sinrux^  I 

F,(')21’1— ! —  w'**0(p) I 
smocx^  I 


aQ  cos  9 


(6.27) 


where  «i,a2’Pl’P2  are  chosen  according  to  (5.9c). 
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Analogically  from  (6.19)  we  have: 


•  F,(-1)  2“'cotn|3,  (s,  *  1)"'  +  F,(1)  2n’cottia.,  (1  -s,)  j 

^■Tf,(1)2"*— w;‘*0(p)l 

2  I  sinna^  I 

1M F,(1)  2"’— !— w!*  ♦  0(p)l 


0 

3(1)2 


p 

F4(1)2 


+  iH^ 


:t  +  1)n’ 

♦  | 

1 

« 

w  : 

sinticx^ 

W1 

1 

« 

w  J 

sinrnx^ 

W2 

1 

« 

w  j 

sinrnx^ 

W3 

1 

a 

W  ‘ 

sinrux^ 

W4 

1 

« 

w 

sintux^ 

W1 

1 

« 

w 

sinna^ 

W2 

1 

a 

w 

sinruxj 

W3 

,  1 

a 

w. 

ctq  cos8  sinB 


(6.28) 


Similar  operations  with  (6.25)  and  (6.26)  when  terms  with  f-j^  are  singular  at  pi=2^=0 


or  x-|=-l  and  S2=1  give  us: 
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33|  "F3(-1 )  2**cotrtP2 (s2  +  l/1  ♦  F3(1)2f’1cotna2(1  -s/* 


tC34[-F4(-1)  2*‘ cottiP2  (s2  ♦  1)°*  ♦  F4(1)  2'  cotrux^  (1  -  s2) 

-F1(-1)2*1— (-2,)*'  ♦  0(p)l 
2  I  1  sinn-p,  1  I 

.  lM  -  F,(-1)  2*1  — 1— (-  z2)Pl  ♦  0(p)l 
2  I  1  sinnp,  *  I 


♦  !5?l  *  F.(*1)  2 


*»  1 

sinnp, 

*1  1 

sirmp, 


(-2,)'  ♦  0(p) 


1,  1 
(■  2„)  ♦  0(p)  J 


.  1^?|  -F2(-1)2*‘— !— (-2,)"'  ♦  0(p)  I 
2  I  sinap,  1  1 

.  (-2/’  *  o(P)  1  = 
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nC43[-F3(-1)  2  1  cotn{52 (s2  +  1  )*  +  F3(1)  20|cotna2(1  -s2)"*J 
+  TlC^[-F4(-1)2*,cotnp2(s2  +  1)°1  ♦  F^(1)  2P*cotna2  (1  -  s2)"*J 


4  ^jrl  -  F,(-D  2”'  — "—(-2,)'  *  °(P) 

0  1  sinrtp,  J 


*TL3[-F'(-1)2',^(-^'*0(p-] 
■  4  0(p)] 

^{•F^1)2“'ii;(-^'*0(p)] 


♦  -F2(-1)2“’— ! T(-z< 

2  I  sinap, 


Ct  1  i 

T-f-z,)'  *  °(P)  =  0 


(6.30) 


Substituting  back  for  W1...W4  in  (6,27)  and  examining  behavior  at  x,j  (ie.  multiplying 
(6.27)  by  lim  (1-s-|)'al  where  -1<ai<0  )  we  have:  lim  (1-s,)  1  -1  <»,  <  0 

s,  -*  1  ‘ 

F2(1)  2'  cotaot,  =  0 
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1 


or 


cosnotj  =  0 


“,=-2 


Multiplying  now  (6.29)  or  (6.30)  by  lim  ($2  ♦  1) 


S2-»-1 


-1  <  p2  <  0  we  obtain: 


-  F3(-1)  2  2  cotn(32  *  0 


or 


cotrip,  =  0 


Let’s  examine  (6.27-6.30)  at  s-|-»-1  and  S2~»1  by  multiplying  (6.27)  and  (6.28)  by: 


lim  (s,  +  1)  and  (6.29-6.30)  by: 


v-0 


lim  (1  -  s2)  where  p1=<x2  =  p 


S2-*  1 


We  obtain: 


1  ...  iH9  _  ...  1  1 


•  «C,2F2(-1)  ^.comp  2  ■jvv^sin)ip 


F,(1)-L  — L-[cose.isine(a*ib)f(--i)  • 


ir’F»(1)i3Sp,co*e*',in8(c*M)f^) 


*^(1)7?^|cos8-isine(a-b)|n(-3 


^5F3(1)7?ilcos0-isine(c'id)f('? 


iH 


10r 


^1)i^|COS0*iSine(a*to)|,(-^ 


2 

i  H 


*  1  ^  7f  simiji *cos 0 "  'sln0(a ' ib)l”  ^  ^ 


i  H 


16 


F,(1)-^  — —  [cose  -  isine(c-  id)f(-4)  =  0 


2  41  '  J2  sirmp 
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F2(-1)(-rcC12  cosnp)  +  F30)(-— )  y  [-  Hg[cos8  +  i  sin0(a  +  ib)f 

-H^cosO  +  i  sin0(c  +  id)f +H13[cos0  -  isin8(a-  ib)  f 
+  H15[cos0  -  i  sin0(c-  id)  ]°  J 

*  %(1)(~)  £  [-  H,0[cosS  .  i  sin6(a  *  ib)f 
-  H12[cos0  +  i  sin6(c  +  id)f 

+  H14[cos8  -  i  sin8(a  -  ib)  ]°+H16{cos8  -  i  sin0(c  -  id)  f  J  =  0 


(6.31) 


Analogically  from  (6.28)  we  obtain: 


F,(-1  )(-tx.C21  cosnp)  +  F3(1)(-— )  y  [-H1?[cos0  +  i  si.n0(a  +  ib)f 

xa  /  L 

-  Hig[cos0  +  i  sin8(c  +  id)f  +  H21(cos8  -  i  sin0(a  -  ib)  ]° 

+  H23[cos0  -  i  sin0(c  -  id)  f  J 

*  F4(1>(t0  1  ['  H"[C0S9  *  1  Sin0(a  *  ib)1" 

-  H2O[cos0  +  i  sin0(c  +  id)f 

+  H22[cos8  -  i  sin8(a  -  ib)  f+H24[cos0  -  i  sin0(c  -  id)  f  J  =  0 


Examination  of  (6.29)  gives  us  the  following: 
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tiC33  F3(1)~c°tnp  «'C34f:4(1)-yjcotllP 

*  "^F,('1)7r^p|cose'isine“'^'^ 

^’F'(-1)7?^,cos0“s,ne“/('y 

-  !rJF'(-1)^iIcos9*isin0“'f(-^ 

-!r-5F’(-1,7?li|cos9tislne“/(^)” 

*1T?F‘('1)^IiNf,c“e'isin9"^'^ 

*T1!F^1)^^lcos9-iSin9“2f('^ 

-ir<F^1,7?^,cos0*isin8“/('y' 

-t^T?  i|cos9tisin9“^-y” 


/  x  \  ^  r 

F3(1)tiC33cosaP  ♦  F,(1)nC„COSnP  *  F,(-1)  i  (-  ji)  [  G,  1  COS0  -  i  Sln8  w,J 

d 

+  Gn[  cosQ  -  isinQ  w2f  -  G13[  cosB  +i sin0  w,f  -  G,5(  cose  +  isinBw/] 

/  x  \  ^  r 

+  F  (-1)i(-— )  G10(  cosS  -  i  sin0  o),f  +  G)2[cos6  -  i  sine  <■>/ 

2  L 

-  Gh[ cos8+  i sine co  1  jn  -  G16[cos0  +  isin6u)2f  j  =0  (6 


(6.33) 


and  analogically  from  (6.30)  we  have  the  last  equation: 


A  60 


NADC-871 33-60 


F3(1)nC43c°Sfip  +  F4(1)  aC44cosfip  +  Fi(_1)y(’^)  [gi?(  cos0  -  i  sin0  co/ 

d 

+  G,g[  cos0  -  i  sin8  co2f  -  G21[cos0  +  isin0co1f  -  G23[  cos8  + i  sin8  ai2f 
+  (yO  [  G18[cos0-  isinSwif  +  G20[  cos0  -  i  sin0  w2f 

d 

-  G22[  cos0  +  i  sin0  u)1]p  -  G24[cos0  +  isin0a>2]pJ  =  0  (6.34) 

The  set  of  equations  (6.31-6.34)  gives  a  system  of  four  equations  with  four  unknowns 
F-|(-1),  fr2(~1  )■  fr3(1)-  F4O)  and  the  right  hand  side  equal  to  zero.  We  can  write  the 
above  set  in  following  form: 


V  *  •  -- 


where: 
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2 

-  fl  COS  tip  C12C2jC33C44  -  —COS  tip  ^2^34^23^41 

2  2 
—  cos  tip  C12C43A24A31  +  -^-cos  tip  C12C33A24A41 

2 

+  -q~  COS  tip  0^2 ^44^23^3 1  +  ^  COS  flP  ^^2^21  ^34^43 

tl2  2  1 

+  -^-COS  tip  C21C44A 13A32  +  ^0  A13A24A31A42 

1  n2  2 

T6A  13^24^32^41  "  ■4-C0S  ^21^34^13^42 

2 

tl  2  1 

-  -^j-COS  tip  C21C43AhA32  -  J0  AHA23A31A42 

1  n  2 

+  ?6AHA23A32A41  +  "4"008  rtP  C21^33A14A42  =  0 


a13=  -  H9[cos0  +  i  sin0(a  +  ib)f -Hn(cos0  +  isin0(c  +  id)f  ♦ 
+  HJ3[cos0  -  isin0(a-  ib)  ]°+  H15[cos0  -  i  sin0(c  -  id)  ]n 

ah  =  '  H1O(COS0  +  '  sin0(a  ♦  *>)f  -  HJ2[cos0  +  t  sin0(c  +  id)f  + 
+  Hh[cos0  -  i  sin0(a  -  ib)  f+H16[cos0  -  i  sin0(c  -  id)  f 

a^3  =  -H,7[cos0  +  i  sin0(a  +  ib)f  -H19[cos0  +  i  sin0(c  +  id)]0  + 

+  H21[cos0  -  i  sin8(a  -  ib)  f +H23[cos0  -  isin0(c  -  id)  ]n 
a^4  =  -  Hlg[cos0  +  i  sin0(a  ♦  ib)]°  - H2O[cos0  +  i  sin0(c  +  id)]0  + 
+  H22[cos0  -  i  sin0(a  -  ib)]n  +  H2<,[cos0  -  isin0(c  -  id)  ]° 

A31  =  G9[cos0- isinOco^0  +  Gn[  cos0  -  i  sir>0  u)2f  - 

-  Gn(  cos0  +i  sin0  u^f  -  G15[  cos0  + i  sin0  w2f 


(6.36) 
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A32"  G10[  cos8  *  i  sin9  cd,]0  ♦  G12[cos0- isin8u)2f  - 

-  GH[  cos8  + i  sinS  -  G16[cos0  +  isin0o)2]( 

A41  =  G17[cos0- isin0wtf  +  Gw[cos0- isin8w2f  - 

-  G21[cos8  +  isin0  w1f  -  G23[cos8  +  isin8cj2]f 
A42=  G18[cos0- isin0o)1f  +  G20[  cos8  -  i  sin0  w2]p  - 

-  G^cose  +  isineu^f  -  Q2A[  cos0  +  isin8a)2f 

The  only  unknown  in  (6.36)  is  P  which  has  to  be  found  so  (6.36)  is  satisfied.  Once  p  is 

found,  two  equations  from  (6.35)  can  be  chosen,  that  together  with  (6.11-6.12)  become 
the  additional  conditions  necessary  to  solve  system  of  integral  equations 
(3.31  -3.34)  by  Gauss-Jacobi  or  Lobatto-Jacobi  method. 

It  should  be  pointed  out  that  in  order  to  avoid  multivaluedness  of  solution  only  one 
branch  of  complex  number  must  be  taken  for  the  calculations. 

The  method  of  choosing  the  branch  is  explained  as  follows.  Formula  (5.9)  is  true  for 
representing  crack  a-b  and  complex  number  z-a^e'^1  and  z-b=r2ei(!)2  where 
argument  is  determined  uniquely  by  making  a  cut  of  the  complex  plane  between  poles 
a  and  b  so  that  0«J>i  <2ir  and  -tc«}>2<tc.  In  this  way,  any  complex  number  z  taken 

to  real  power  p  will  have  a  unique  solution  as  follows: 

[cos8  -  isin8w,f  =  pg  [e* ]p  -  pg  e* ^ 


'/  V  .-'V.- V  ’ji  - .  v.^.. 


■:v,v/Wav.v.'/.;-! 
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VII.  Stress  intensity  factors  for  skew  crack  case. 


The  set  of  singular  integral  equations  with  the  generalized  Cauchy  kernels  and  four 
additional  conditions  in  the  form  of  the  equations  (6.11),  (6.12),  (6.31)  and  (6.32)  can 
be  solved  by  Gauss-Jacobi  collocation  method.  We  can  find  t^  as  roots  of  Jacobi 


polynomial  of  order  n  as  follows: 


p'.T'V)  -  o 


Collocation  points  Xj  are  found  as  roots  of: 


x  j  .  0 


and  weighting  coefficients  are  calculated  according  to  formula: 


Wk  =  * 


_ (2n+ctn>j3+2)J[Xn+a+1^ _ 

(n+1)-  (n+ot+p+1)  r(n+oc+p+1)  p£f(  g  |p^p)(  g 


In  order  to  find  expressions  for  the  stress  intensity  factors  we  start  from  definitions  (5.1- 
5.8)  and  follow  the  same  procedure  with  different  order  of  singularity. 

Starting  from  (5.10)  with  X5=0. 


ki(xJ  =  lim  +s2  ay,y,(s2'°) 


Since  at  xa  only  the  two  first  terms  of  (3.33)  are  singular,  using  (5.9)  we  have: 
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m  **  k  %“ 


Substituting  the  above  results  for  stress  in  (7.1)  and  after  some  algebra  we  get: 


k,(x.)  -  *  C«F<('1)] 

k2(xj  -  ti^yoT  [C.a'SH)  *  V/'1)] 

k,(xd)  c12f2(1) 

^2(x,j)  = '  n  2B  JX4  ^2)  Ft(1) 


where  rest  of  stress  intensity  factors  were  found  in  similar  way. 
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VIII.  Normalized  stress  intensity  factors. 

In  order  to  present  computational  results  for  two  crack  configuration  and  for  skew  crack 
configuration  ,  the  stress  intensity  factors  are  normalized. 

Let's  define  length  of  the  cracks  as  xaxjj=L-j  and  Xc*d=,-2-  Then  normalized  stress 
intensity  factors  become: 

P=1>2  (8.1) 

Thus,  the  mode  one  stress  intensity  factor,  for  the  two  crack  configuration,  at  xa  is: 

*  C«F«(-1>]  (8.2) 

CT0 

and  for  the  skew  crack  case: 

klOg  =  —  [  C33F3(_1>  +  CMF^’1)]2  (8.3) 

ao 

IX.  Strain  Energy  Release  Rate  for  x^. 

Let’s  use  the  usual  definition  of  the  strain  energy  release  rate  [3]: 

G-^fU-V)  (9.1) 

where  G  is  the  total  energy  and  consists  of  two  parts  which  represent  mode-1  (  G|  ) 
and  mode-ll  (G|| ): 

G  =  G,  +  G(|  (9.2) 
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at  x  =  xd,  we  may  write  [4]: 


dU-dV  _  1  f  f  cr 

da  ’  2  J  1  yiyi 


ffv<y(xi'°^Y(xi‘da'0  )  ’  y(xrda,0')] 


+  Vy  (xi'°)l  u(xf  da'°  )  -  u(xf  da,0*)]  dx, 


-0')l]< 


G,  =  Ija^Ot^IvOCt-da^  -  v(x,  -  da,0")J  dxt 

Ni 

x^  +  da 

Gn  -  ilTx,y1(xi,0)[u(><,'da'0+)  *  U(xrda^°’)]dxi 


To  evaluate  integral  (9.3),  the  asymptotic  expressions  of  stresses  and  displacements 
around  the  crack  tip  are  needed.  The  expression  of  normal  and  shear  stress  can  be 
easily  obtained  by  using  the  following  definitions: 


k.(  x  ) 

c  (x.,0)  =  —  --  — r  +  higher  order  terms 

v,v’ 


v,(xi-0)  ’ 


'2(k,  -  xV 


higher  order  terms 


Asymptotic  expression  for  the  displacements  can  be  obtained  using: 


Au  =  u(xt,0+)  -  u(x,,0’)  «  Jim  ^  J[  cj^  -  c^Jdx, 
AY  =  v(x,,0+)  -  y^x^O")  =  Jirn  ^  J[c^  -  c^Jdy, 
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where  total  strains  can  be  found  from  (1.22-1.23)  from  total  stresses  obtained 
according  to  (3.4a,b)  using  (1.14-1.17)  and  (3. 2-3.3)  upon  substitution  of  (3.1)  into 
(2.20-2.25).  However,  it  is  known  that  at  x<j  only  terms  with  fi  and  I2  contribute  to 
singular  behavior.  Thus  Au  and  Av  can  also  be  found  directly  us.  .g  (1 .20-1 .21 ),  [5]. 


AU  =  |f1(x,)dx| 

(9-8) 

AV  «=  jyxjjdx, 

(99) 

Starting  from  this  point  we  shall  substitute  for  f-j,  I2  from  (1.29  -  1.30)  taking  A,B  from 
(1.35-1.36)  as: 

A  =  K|  t 


B  *  -K1  -  K2tOj 


where 


K,-l - 1 

1  i 


2^,2  2.  J 

s  2an(u>1  -U)2)  ^ 


I1 w 


e  'dt, 


K2  " 


|s|  W1W2 

.3,2  2 .  j  v  ^ e  dti 

I  s  2a22(w1  -  «2 )  ^ 


Hence,  the  normal  displacement  expression  is: 


+*•  2  -«Js|y.*  «i|s|y] 

-1  r  s  Ae  +(Ac.+  Bc2)e  2 

ay  =  lim  —  rr  _ 1 - 1 - ra  </ .  a 

V)-o  2rtilsl  w  ^  12  1  22; 


„  -«,|s|V:  ^  «2|s|yJ 

Be  2  1  ♦  (Ac.  -  Be.  e  2  , 

*  - 2 - ! -  (»12«  -»22) 

“2 


-ISX. 

e  n  ds 


(9.10) 
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and  for  shear  displacement  the  analogical  expression  is: 


Au  =  lim  [  -• 

y,  ->  o  1 


(Ae  -  (Ac ,  +  Bc2) e  ](anu>*- a,2) 


r_  *®2ls|yi  «Js|y;  , 

4  (Be  -  (Ac3-Bc,)e  l(a„Uj-a,2) 


*'SXt 

e  nds 


(9.11) 


where  c-\  ,02,03  are  taken  from  (1 .13  a,b,c). 
With  (9.12-9.15)  we  obtain: 


Ay 


=  lim  ~\ 


1 


v 2|s|(^-^) 


\w  *S\ 


'11 


2  . .  -«,WW  «,|s|yj 

(at2wt  -  a22)(e  -  e  ) 


2  ...  "^ivi 

(at2C02  ‘  a22>(e  ) 


+  Jii^L 


Jyv6'*'' 


dt, 


IS  a22  |(at2wt'a22)b>2 


-®,|s|y(  «,|s|yj  2  2 

(e  -  e  )  -  (at2w2  •  a22)u,((e 


'u2)sK  “zislvj 
+  e  ) 

. 


) 


-ISX, 

e  ds 


(9,12) 


and 
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For  skew  crack  configuration 

xc  =  °  t,  =  l~xd  dt,  =  lxddT  -1<t<1 
the  equations  for  displacement  now  become: 


Ay  =  lim 


i  If?*, 

e  2  dr  ffm2  -  —  V 


ls|yi, 

P  1 
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-  ~  f  *  ^ 

..  *xd  f  *  2  f, /N  2  \  u  2  ai2w‘"iWtf  Atti|s|yK 
Au  =  hm  - r — r  ft(T)e  dT[(a)r  — )(e  +e  )- 


v,->°  -  is(“i  *w2)  -i 


1  IfV 


*  *e“  ’)]  4  — — -  ff2(T)e  2  dr[(0)2--^)w 

ai1  ISW1W2-1  aH 


^-“1|S|W  “llslyK  /  2  ai2,  ,  -"zHW  “aMvL  -isxl  .  /Qi.. 

(e  -e  )  ’  (U)2T')U)i(e  'e  )]  e  ds  (915) 


First  we  have  to  evaluate  integrals  of  terms  with  variable  x.  Note  that  they  all  have 
general  form: 


■IVD.-* 


fe'at(1  +T)pFj(T)dr  (9-16) 


where  a  = 


(9.17) 


It  was  assumed  that  function  Fj  can  be  represented  by  the  following  series  expansion: 


Fi«  -  2  cl!' 


(9.18) 


Now  after  substituting  this  for  Fj  a  series  of  integrals  is  obtained 


-  Vo -ro  -tjndt 


Applying  the  Rodrigues'  formula  and  integrating  by  parts  it  reduces  to: 


_  (-1)n  f/-  v" 

I2  =  -  (i  a)  e  (1-t)  (1+t)  dr 

(■2)  n\i 
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Now  upon  transformation  t  =  2u  -  1 ,  dx  =  2du,  0  <  u  <  1  and  simple  algebra  the 
integral  reduces  to: 

. .  .  n  J 

I2  =  2 e'1*  ei2*“(1  -  u)*“  u"“  du 

2nn!  o 

and  using  the  result  from  [6],  it  gives: 


T  (ia)  -a+o+cn+i  -i&  r(n+(3+1)  r(n+ot+1)  .  _  „  .  ^  .. 

I7  =  - — —2  e  — — 1 — $(p+n+1;2n+cx.+p+2;  2ai) 

2nn!  r(2n+<x+p+2) 


(9.19) 


where  <t>(a;c;z)  =  iF-|{a;c;z)  is  a  confluent  hypergeometric  function  and 


r(x)r(y)/r(X+y)=B(x,y)  is  a  Beta  function.  The  series  expansion  is  true  for  any  x  from 
interval  (-1 ,1),  so 


n=1 


(t)  F,(1)  Ft(1)  F,(1)n!r(a+1) 

c  = - m - n!  = - 

n  ,(«,!»), 


P^a'n)(1)  (a+1)n  T(cx+1+n) 


since 


a=T 


and  analogically 


n!  r(y) 

s^-FiO) - r-  where  r(y)=2v^ 

r(n  +  i) 


n=1 


(2)  F20)  f20)  ,  F2(l)n!r(ofl) 

c  = -  =  - n!  = - 

n  ,(«,!>), 


pKP)m  (a+1)n 


r(<x+1  +n) 


(9.20) 
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Applying  all  these  results  to  equations  for  Av  and  Au  yields: 


f  s" 

av  =  lim  A  1-7  $(p+n+1;  2n+ot+p+2;  isx,) 

yi  L  Is' 


(I)r/  2  ai2  .  ,  -“aW  «^iv; 


Cn  '!K~H(e 


) 


‘11 
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By  changing  the  integration  interval  to  (0,»),  the  above  expressions  for  crack  surface 
displacements  become: 
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We  now  apply  the  asymptotic  expansion  tor  confluent  hypergeometric 


functions: 


_  .  r(c)  e™  r(c)  x  *-c 

$(ajc;x)  =  — —t—  Z —  x 

r(c-a)  _  x  _  r(a) 


(9.24) 


c  =  1  if  lm(x)  >  0 


c  =  -1  if  lm(x)  <  0 


where  in  our  case  a  =  p  +  n  +  1;  c  =  2n +a  +  p +2;  x  =  ±isx,; 


Finally  applying  this  to  the  equations  for  displacement  we  obtain: 
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Foliowing  the  integral  3.381  from  "Table  of  integral"  by  I.S.Gradshteyn  and  I.M.Ryzhik 

[71: 
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we  obtain: 


°  +  2 


ay  =  -  2  F,(1)yxd(xd-x)'  +  higher  order  te .ms 


and 
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for  xi  <  x,-j 

Now  ,  substituting  in  (9.4-9.5)  and  integrating  we  have: 
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or  substituting  for  C-\2  and  C21  and  computing  the  total  strain  energy  release  rate  we 
have: 
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(9.32) 


which  has  the  form  reported  by  Cherepanov  for  the  orthotropic  case. 

For  isotropic  material  and  self-similar  crack  propagation  this  reduces  to  the  well  known 


form: 


G  = 


nk^(xd) 


(9.33) 
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